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It is frequently convenient in data processing to convert analogue 
to digital data for computer assimilation. A convenient method of such 
conversion has been developed and used in the study of correlation detec- 
tion of audio signals corrupted by noise. 

A method to use apriori knowledge of the corrupting noise to increase 
processing gain has been studied. In the case of detection of a sinusoid in 
noise, an additional gain over conventional auto-correlation of up to 14. 5 
db has been achieved. 

Finally, a signal source located in an unknown random noise field 
was detected, classified and located in relative bearing by the cross-cor- 
relation of the signals received from two spatially separated sensors. 
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Introduction. 



During the summer of 1965, the author visited the U. S. Naval Elecri 
tronics Laboratory, San Diego, for six weeks. The.. experience gained there 
suggested several problem areas which would profit from further study. 

It appeared that there was a need to analyze large quantities of ana- 
logue data collected during various experiments. The volume of the data 
almost precluded "hand" analysis and y^ no facilities existed to enable a 
computer to assimilate the data. This facility was also lacking at the U. S. N. 
Postgraduate School Monterey. 

The data collected during these experiments at San Diego was often 
corrupted by unavoidable locally generated noise. It was felt that if more 
could be learned about this self-noise some better method of data extraction 
might be possible. 

One of the most powerful methods of data analysis is the correlation 
technique which may reveal information masked by corrupting noise. These 
techniques have been thoroughly investigated over the past 25 years, and a 
wealth of literature exists on this subject. 

Although continuous time integrations and correlations are possible 
mathematically on known explicit functions, physical data often does not 
fit a known model. In this case approximations to infinite integrals etc. , 
must be made using either analogue or digital equipment. Analogue equip- 
ment tends to be limited in versatility and is expensive, while if a general 
purpose digital computer is available, it can be made to simulate the physical 
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situation. Therefore attention was directed toward describing the analogue 
data by a digital computer. The use of a digital computer to compute various 
correlation functions is the subject of this report. 

2. Statement of the Problem. 

The correlation function of two continuous time variables is given by 



where s(t), N(t) are defined for -«©<t< « 

When the functional form of s(t) and N(t) are known it may be possible 
to evaluate Sq. (l)mathematically. When the form of s(t) and N(t) are known 
only as analogue signals Eq.(l) may be instrumented by a multiplier and an 
integrator. 

The same effect may be achieved by digital means, subject to errors to 
be discussed. 

The sampling theorem states that a continuous function can be repre- 
sented in a finite interval by a finite number of samples of the function. The 
rate of sampling must exceed twice the highest frequency component of the 
Fourier transform (spectrum) of the function. 

By suitably band-limiting a function, samples may be taken at a prac- 
tical rate. The samples if digitally described can then be processed. Equa- 
tion (l) implies a memory extending into the infinite past which cannot be 
achieved physically. As an approximation the "short-time" correlation may 
be formed with some resulting error. 




( 1 ) 
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If the "short-time" correlation is formed digitally, errors arise from 
both the limited number of samples and by the finite accuracy of the indivi- 
dual samples. If "short-time" is defined to be some period much longer than 
the period of the lowest frequency in the signal, errors of the first kind above 
tend to zero as the length of record grows. Measurement error is negligible c; 
if the samples are described to four place accuracy. 

The "short-time" correlation is formed as 

N-K 

=1 < 2 > 

i = 1 

where 



AT - the sampling interval 

, N?' = sample values of the continuous functiontaken at 
intervals of T 

k = o, 1, 2 K 

N = the number of samples available. 

In the limit as N goes toe>o , andA T 

^S*+ 

For finite N 

A 

Rg*+j^*(kAT) = (4) 

Equation (2) holds only if the data has zero mean. If in Eq. (1) we 
assume that s(t) is the desired signal, of zero mean and that N(t) is random 
Gaussian noise and that the observed V(t) = s(t) + N(t), then we may form 
T). This function may in turn be broken down into components: 
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Since the noise and signal are assumed independent the third number 
vanishes in the limit T . 

We may approximate Eq. (5) as 

in which we wish to determine Rg*g*(k^T). 

For very large k = Kq the second member of Eq. (6) will tend toward 
zero because of independence leaving 

%*+N*^^o “ Rs*S*^^o ^ 
j = o, 1, 2, J 

4 ; 

This restriction is not desirablea since if Rj^*j^*(kAt) were known for 
all k, then in Eq. (6) 

A. A. A 

= Rs*+n*( 1^^T) - Rj^*j^*(kAT) (8) 

However Rj^*j^*(kAT) is never known exactly. It can only be known, 
since it is a random variable, in an average sense. 

If we assume that the normalized correlation function Rj^j^(kAT) is 
known apriori 

Rnn(I^^T) = RJ/|^J(kdT) (9) 

where is formed in the absence of signal 

is to be determined where 0 ^ 1 

then in Eq. (8) we may write 

Rs*S*^^^^ ~ Rs*+n*^1^'^T) - ^ Rj^j^(kAT) 

Methods of determining have been studied for the very limited case 
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of a single sinusoid signal masked by noise. 

To gain some experience with the technique of cross-correlation, a 
low level signal was placed in a noise field. The putput of two spatially 
separated sensors was processed to detect and locate the signal. The re- 
sults are presented in section 5. 

3. Conversion of Analogue Signals to Digital Form for Computer Analysis . 

Maa^ analyst of real data may be performed on a digital computer by 
first manually sampling the data a sufficient number of times, and entering 
the sampled data into the computer. 

When .the number of data samples desired exceeds a few thousand, some 
automatic means of data conversion must be used. 

The method should preferably be "off-line, " that is, use auxiliary 
equipment, for efficient employment of the main high speed computer. The 
Computer Laboratory of the Electrical Engineering Department was equipped 
with the necessary items of hardware and these facilities were used to sample 
the analog data and form the digital record. 

The three main items used were a 12 bit analogue to digital convertor 
a Control Data Corporation (C. D. C.) Model 160 Computer having a capacity 
of 4096, 12 bit words, and a C. D.C. Model 163 magnetic tape unit. 

As now constituted the analogue data is consecutively sampled by the 
analogue to digital converter working under control of the CDC 160 computer. 

When 4000 samples have been taken, the computer writes the samples 
on magnetic tape in a format that is recognized by the main CDC 1604 comp 
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puter. After a delay of 0. 4 second another 4000 samples may be collected, 
and so on. Each group of 4000 samples when written on tape is accompanied 
by an identifying word so that any block of 4000 samples may later be located. 

The maximum sampling rate is 5 Kcs. implying that all signals shopld be 
band-limited to less than 2500 cps. The sampling rate may be varied at will 
by manual program entries, as can other convenient parameters. 

The system is capable of digitizing only one source at a time. To allow 
cross -correlation of two or more sources, a clock pulse input is used to in- 
itiate a round of sampling. 

At present a simple AND gate is used to sense the presence of a clock 
pulse. A more permanent arrangement should make provision for the detection 
of a simple clock code, to eliminate the occasional false starts due to noise 
pulses, which have occured under present arrangements. 

Once the data has been recorded digitally on tape, the data may be 
read into the main computer using standard programming techniques. The 
subroutines to accomplish this are available as Subroutine Data from the 
Computer Facility. 

The data when read into the computer is delivered in integer format, 
i. e. , without a decimal. To convert it to units of volts each data point is 
then divided by 409. 6, a factor inherent in the particular A/D converter used. 
The diagram below shows the essential steps. 
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Figure 1. Flow Diagram of Digitizing Process. 

The system is presently adequate to describe many physical processes. 
Certain signals, such as speech, are beyond the capabilities of the system. 
However modifications to the Computer Laboratory hardware to improve per- 
formance were not attempted because of pending contractual arrangements to 
completely replace the hardware with more sophisticated equipment. The new 
equipment will enable two records to be digitized simultaneously. The memory 
capacity will be increased by nearly a factor of four, and the sampling rate 
will be approximately 20 Kcs. This equipment will^be available for use in 
mid 1967. 
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The author has recorded the mixed output of a signal and a noise gen- 
erator to provide a vehicle for the study of correlation. Several other students 
have used the system for their own particular needs with no difficulty 
experienced to date. Additional details on the digiting process are presented 
in Appendix A. 

4. Noise Removal in Auto -Correlation . 

Correlation techniques have gained considerable favour in recent years 
as a method of signal enhancement. Theoretical derivations of the processing 
gain to be expected have been stated by Lee[[l J, for both auto and cross- 
correlation. 

Auto and cross -correlation are defined as follows: 



where 



where 



R=(r) = lim 






T/2 



v(t). v(t + r)dt 






-T/2 



V(t) = s(t) + N(t) 

s(t) = the signal to be recovered 

N(t) = noisaj. here assumed Gaussian 
1 1 

Rc('r) = lim V(t). sMt + r)dt 

T-^ ^ -T/2 



( 1 ) 



( 2 ) 



s^t) is chosen apriori to have the form of s(t) and is locally 
generated . 



Approximations to Eq. (l) and (2) may be generated using sampled 
data in which case. 

N 






(3) 



i = 1 
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( 4 ) 



(5) 



) = W ^ Vi(i) - V2(i + r ) Rc(T) 
i = 1 

A 

where Vi = V(t) for t = iA.t 

It may be shown that Eq. (3) and (4) are of the form 
= Rss(T) + W'l”) + 

Of the four terms in Eq. (5) it is desired to determine only the first, 

A 

Rss(T) which represents the signal contribution. The remaining terms are 
undersirable and success is measured by the extent to which they can be 
eliminated. 

A A 

The cross-terms Rgj^, Rj^g will vanish as N— ><*>, because of the as- 
sumed independence of the signal and the noise. 

Since the correlation is to be done using a finite number N of samples, 
these terms will become small, but may be significant. 

Let us assume that a signal is of the form 
S(t) = TT E cos wt 

N(t) = Gaussian noise of zero mean and variance 6'^ 

The variance of the output of Eq. (3) and (4) has been evaluated by 
Lee LlJ, who clearly shows the superiority of cross -correlation over auto 
correlation. 



In the case of auto -correlation 

2 ^ N ^ N 



while for cross -correlation 
*2 _ 



% = 



+ E^CV 



2,r2 



N 



( 6 ) 



(7) 
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This improvement is achieved only when the form of the desired signal 
is known apriori. In many cases of interest particularly in the passive sonar 
situation^ the desired signal form is not known apriori, and consequently 
auto-correlation must be used. 

In addition to the noise input supplied by the environment, local dis- 
turbances may cause noise effects to appear as part of the correlator output. 
If these two sources of noise, taken together, are known apriori, and can be 
considered to remain stationary over the interval under consideration, then 
perhaps some means can be found to eliminate their corrupting effect on the 
output of the auto -correlator. 

This problem has been studied for many years. In an early work 
Lee states 



". . .the correlation of the noise, once measured 
can be compensated, (for)'' 

The details of the implementation of this were not specified. 

In a much more recent work Van Trees [sj considers an optimum 
receiver in which he considers that the noise is known apriori as a norm- 
alized correlation Y (T). After correlation of the received data some 
constant K times / ( 'V) is substracted form the correlator output. Deter- 
mination of K rests upon knowledge of signal strength and noise strength. 
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Again it is not clear how to instrument the calculation of K. 

We will assume that the sample auto-correlation of the signal. 



Rs+n(^) has been formed in the presence of corrupting noise. We shall also 
assume that the statistics of the corrupting noise are known apriori and are 
represented by the normalized correlation function formed in the 

absence of signal. 

The object of the processing is to determine a linear combination of 

A ^ 

which will enhance the signal and tend to remove the 
noise. No claims are made for the optimality of the method presented here, 
other than it does improve the signal to noise ratio of the output. 

In the following we will be dealing with sampled data functionS/not 
continuous functions of time. To simplify notation however the symbol T 
will be used to denote a sample data function. 

The power spectral density P(w), and the auto-correlation functions 
are transform pairs, that is 



Given one member of the pair, the other may be found. Therefore the 
processor could attempt to operate on either the power spectral density or 
the correlation function. 

From Eq. (5) we get (neglecting cross terms). 



R(T) 




( 8 ) 




(9) 
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( 10 ) 



or equivalently by Eq. (9) 

Pgs(w) = Ps+N^w) - 

We have assumed that as apriori information estimates of the norm- 

* 

alized power spectral density function P^^Cw) and normalized auto-correla- 
tion function » formed in the absence of signal^ are available. At 

any given instant there is no assurance that the apriori known noise estimate 
represents the true noise in other than an expected value sense. 

Then in Eq. (lO) and (11) we may write 

RSS(T) = RS+n(T) - Rnn(^) (12) 

Pss(w) = Ps+n(w) - Pnn^'^^ ^13) 

where are known functions 

O $ C I / o< ^ constant to be determined. 

Determination of the ^ in Eq. (12) and (13) becomes the central 
problem to be considered. 

The following methods of calculating ©< hold only for the case of sin- 
usoidal signals in noise. Further study is necessary before any general 
method can be stated. 

Since the correlation function was produced from the data directly and 
the power spectral density indirectly from the data via the correlation, it was 
thought advisable to operate first on Eq. (12) rather than Eq. (13). 

In the case of a sinusoidal signal, it was thought that finding an 



24 



Max Rss(T ) 



(14) 



^ such that Rss(0^ 



would be advisable. 

However Eq. (14) is only a necessary not a sufficient condition for 
the formation of a valid, physically realizable auto-correlation function. 

Thus an instrumentation of Eq. (12) to satisfy Eq. (14) may result in 
a correlation function which is the output of a physically unrealizable linear 
system. In the time domain there seem to be no simple tests for realizabi- 
lity. 

Several signals were analyzed using Eq. (12) and (14) at input signal 
to noise ratios of 0 and -10 db. No trouble was experienced with the 0 db 
case; the signal was almost completely recovered and the noise suppressed. 
The -10 db input signal after similar processing was transformed inl;o the 
frequency domain. The signal was clearly evidenced by a large spike in the 
power spectral density at the signal frequency. The amplitude of the signal 
spike after processing showed nearly 6 db greater discrimination between 
signal and noise peaks than it did before processing. 

The above gain would have been encouraging had not the power spect- 
ral density,' showed some negative contributions. 

Negative contributions in a power spectral density, are evidence of, and 
in fact are a sufficient condition for, physical unrealizability. 

Because of this difficulty no further efforts were made to suitably deter- 
mine on the basis of the correlation functions. Efforts were then directed 

to a determination via Eq. (13). 
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If the signal is a pure sinusoid with no noise present, the power 
spectral density shows a delta -function cT (w) at w = ws, and is identically 
zero for all other values of w. The perfect processor would give this output 
even in the presence of noise. 

There are two conditions one may then place on Eq, (13). These are 

(a) There must be no negative contributions 0“=^ wi: w 

Max 

(b) At least one point must be zero. 

This allows an to be determined directly in the frequency domain. 

c 

Condition ,;(a) above ensures realizability, while (b) ensures that for any 
, condition (a) is violated. On the other hand choice of an 
results in the removal of less noise than one might achieve other- 
wise, but at least ensures physical realizability. 

A trial and error search routine was carried out to find an o< which 
satisfies conditions (a) and (b) above. The search routine is halted when 
ano< is found which results in at least one point of Eq. (13)having a value 

^ .00001, and no negative points. The time to complete a search 

is in the order of . 5 seconds, which is small compared to the time taken to 
form the functions originally (70 seconds). 

Having formed results were renormalized to make 

the area under the curve unity, and then plotted. 

/ 

The determined by the search was then used in Eq. (12) to form 
c 

Rss(T^) and the result plotted. No claims for validity are made for this step 
except that having formed (T) in this mannenits transform had no negative 
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components. 



The noise removal proceedure based on an ck determined from the 
power spectral density resulted in an additional 14. 5 db processing gain. 

This figure was obtained by considering the ratio of energy under the 50 cps. 
comb-width to the total energy under all other comb-widths, both before and 
after the noise removal. Before noise removal the calculated input signal- 
to-noise ratio was -8.84 db. , After noise removal the equivalent calculation 
yielded +5.62 db. 

Another measure of effectiveness should be mentioned. In many cases 
the ability to distinguish a signal from noise in the spectral density is of 
interest. In the power spectrum the ratio of the amplitude due to signal to 
that due to a noise "spike", may be expressed in db. In the -8. 84 db 
(noninally -10 db) case, the ratio of signal to highest noise peak was 5.62 
db. before noise removal, and 9. 05 db. after noise removal, or an improve- 
ment of 3 . 4 db. 

Details of the processing are included as Appendix B, 

5. Location of a Source by Cross -Correlation . 

The bearing of a signal with respect to an abserver may be determined 
by measuring the difference in arrival time of the signal at two or more sep- 
arated observation stations. 

As an experiment a 200 cps. sinusoidal signal generator was set in 
one corner of a cafeteria. Two microphones spaced a distance d units apart 
were located symmetrically with respect to the center line of the room, at 
the far wall of the cafeteria. Interfering noise was supplied by the noon-day 
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clientele. 

The maximum level of signal was limited by the intolerance of the din- 
ers. This was determined to be 38 milli-watts measured at the loudspeaker 
voice-coil. At this level the signal was inaudible at the microphones. 

A stereo recording was made of the signal plus noise, and after slow- 
ing the tape by a factor of four, selected portions of data were digitized 
using techniques detailed in Appendix A. The resulting digital data was 
cross -correlated in the CDC 1604 computer and the correlation function plot- 
ted. 



The geometry of the problem is shown below. 



soufce 



Nfoise field 




70. 7S 



Assuming plane wave propagation, and in the absence of noise, we 
may associate with A D in Figure 1, an equivalent amount of time; that is. 
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A T = A D/C where C is the velocity of sound in air. 

If the two sensors are cross -correlated a maximum in the function will 
occur whenever the time variable of Correlation is equal to a multiple of the 
period, of. the signal. 

The signal was located successfully to an accuracy of + 0. 17°. The 
frequency of the observed correlation function was nearly 150 cps. , or three 
times the expected value of 50 cps. This would Indicate that detection was 
accomplished on the third harmonic of the source oscillator. 

A study of the noise characteristics supplied by the diners was not 

made. 

Provided that accurate time synchronization between the two data 
sources to be correlated can be maintained during the digitizing process, 
correlation will provide both bearing and frequency information. 

Details of the problem are presented as Appendix C. 

6. Conclusion. 

A convenient method of digitizing analogue data has been developed 
and this digital data used to detect signals masked by noise. 

The digitizing system is now operational. 

Many interesting areas of research remain to be pursued in the field 
of recovery of noisy signals. In particular additional effort is needed in the 
area of removal of aprlori known noise effects. 

The facilities to perform the necessary calculation on a digital computer 
are now available at the U. S. Naval Postgraduate School. 
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APPENDIX A. 



ANALOGUE TO DIGITAL DATA CONVERSION 

1. INTRODUCTION. 

The process of translating real, analogue data into computer digitized 
form has been accomplished in two separate stages. 

In the first stage, the analogue data on magnetic tape is sampled at 
discrete points of time. These samples are digitized and written in bit form 
on magnetic tape. 

In the second stage, the digital magnetic tape is loaded on the CDC 
1604 peripheral tape units and then read into the computer. 

Earlier work at USNPG School in computer analysis of analogue data 
utilized the "on line" approach, in which the data was fed to the main com- 
puter directly from the digitizing computer. This mode of operation is now 
impractical because of high computer facility utilization. 

The method detailed here results in data being reduced to a form suit- 
able for submission as a standard program input. 

2. DIGITIZING THE DATA. 

Basically the method takes samples until the Core storage of the dig- 
itizing computer is full, and then writes this information on magnetic tape 
in digital form. 

The digitizing computer is a CDC 160 located in the Computer Labor- 
otory of the Electrical Engineering Department. 

This computer has a 4096, 12 bit memory length, and has the ability 
to select various external equipments through the proper external function 
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codes. The lower 4000 cells of the computer were reserved for data storage 
and the top 96 cells for programming. Thus the core storage available sets 
one quite restrictive limitation on the system. In practice however, 4000 
samples were adequate to describe many common processes. 

The block diagram of the system is shown in Figure 1. The 160 selects 
either the A/D converter, for data input, or the 163 tape unit for data output. 
When the A/D converter is selected (DCF 14xx) and the INA (IN put to A) in- 
struction is executed, the input analogue signal, is sampled, and a 12 bit 
number representing amplitude is sent to the 160 and stored. Repetition of 
the process in quick succession allows blocks of up to 4000 samples to be 
stored. The maximum sampling rate achievable is 5000 samples/second. 

To control exactly the time at which samples are taken, the execution 
of the INA instruction is prevented until two events have occurred; first, an 
enabling switch must be closed to ground (EXF -2410), and second, a pulse 
(see Figure 2) must be present at the base of the AND Gate» When these 
conditions have been satisfied the flip-flop goes to the "1” state placing -3 
volts on the INA line and thus executing the INA instruction. 

Further samples may now be taken as the flip-flop will stay in the "1" 
state until EXF 2400 is called, resetting FFI to the "O" state. The complete 
process is illustrated in Figure 3. 

The usual mode of operation is to have each block initiated by a pulse, 
with the intersample delay under computer control. Under certain circumstan- 
ces, it may be desired to have only the runs initiated by a pulse, with both 
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the block count and the sample count under computer control. 



The inter-block timing in this mode is not accurate enough for cross- 
correlation. The basic program shown in Table 1 may be modified using 
Table 2 to give various modes of operation. 

Operating Instructions-Digitize Eq. (3). 

Features of Digitize Eg. (3) . 

1. Manual Entries 

a. 0066 - intersample delay. See Figure 6 

b. 0067 - any identifier you wish 

c. 0070 - Initial run number 

d. 0071 - number of runs desired 

e. 0072 - number of samples/block ^ 7640g = (4000) 

f. 0073 - number of blocks/run 
Note that (d) x (f) ^ 7600 



2. Main Points 

a. A block of n, n ^ 4OOOj^0 samples is taken each time a clock 
pulse is present and EXF 2410 is made. 

b. The inter-sample spacing interval is program controlled and is 
variable from approx. 190 jAsec. minimum to about 50 Tnsec. maxi 
mum. 

c. The first four (4) words written on the 163, of each block, are 
identifying words in the following order: 

(l). Run number. 
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(2) . Spare (may be preset to anything) 

(3) , Block number 

(4) . Number blocks/run 

d. The identifying information when digitizing data must be recorded, 
as 1604 subroutine DATA uses this information to locate blocks of 
data on digital tape. 

e. At the end of each run a check is made to ensure that the storage 
capacity of a full roll of 163 tape is not being exceeded. If it is, 
an End of File is written and the program halts with the total numb- 
er of blocks written in the A register. 

f. Program may be stopped at any point to allow cueing data tape, 
provided that it is not cleared. Progress resumes when 160 placed 
in run condition. 

g. Analogue, input is at A/D Channel 1. A voltage reference is est- 
ablished by an external source at about -5 volts. Correct opera- 
tion is assured if the proper test pattern is displayed. See in- 
structions for Program Test 160. 

h. A voltage of V volts at A/D Channel 1 is converted to a digital 
reading by the A/D converter, equal to -409. 6 x V 

3. Operating Instructions 

a. Load from cell 0000 

b. Run from cell 0000 

c. When digitizing is complete program stops at 0124 with current 
run number in A register. 
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d. Clear and run from cell 0127 to write E. O.F. 



4. Modification 

a. See Table II for various possible modifications. 

The limitations of the sampling process are graphed in Figures 4 and 5. 
Figure 4 shows the effect of sampling speed on the percentage of the data 
that may be digitized. The speed of the data gathering is a variable under 
program control, while the speed of output is not. Thus at slower sampling 
rates, the sampling approaches continuous sampling, and the amount of the 
real time data which may be represented on tape, iitcreaises;, as seen in 
Figure 5. 

Figure 6 shows the relationship between the number entered as inter- 
sample delay timing word, and the resultant inter-sample delay. 

The development of a clock pulse from initial, recording to output at 
the flip-flop is shown in Figure 2. Experiment has shown that a total delay 
of 100 micro-seconds exists between tQ, the time of initiating the pulse on 
record and t^, the time the pulse is delivered. This ^ period remains station- 
ary over a large range of level settings, pulse widths and amplifier gain 
settings o Thus the clock pulse may be used to synchronize two p^ees.-6f data, 
as in cross -correlation. 

Experience has shown that since only the first of a series of adjacent 
clock pulses in used to initiate sampling, the recording of a train of pulses, 
rather than a single pulse is a great aid in cueing the analogue tape. The 
length of the train must be less than the length of time taken to sample one 
block. 
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In order to set the CDC 160 up properly for digitizing, a short program, 
TEST 160, has been devised (see Table UI). In this program, the first cells 
-0000 to 0017, contain a program to accurately set the bias voltage. With 
no signal in and the bias at -5. 0 volts, the value of the conversion should 
be 0000 or 7777. By monitoring D/A Channel 1 with a "scope", deviations 
from the zero condition may be detected and adjusted using the voltage con- 
trol on the bias source. Small errors in zeroing may be expected, due to line 
drift, etc. The effect of this may be eliminated by removing the mean of the 
data in the 1604. 

The cells from 0020 to 0045 contain a program which has identical 
timing to that of DIGITIZE. By observing the input request pulses on the 
patch box, between the CDC 160 and the A/D converter, one may adjust the 
value of the timing control to give the desired inter-sample delay. Chang- 
ing the control word by one number changes the sampling interval by 12. 8 
microseconds. 

In order to distinctly label each block of data written on digital tape, 
an^ identifying series of numbers is first written, followed by the data points. 
Since it takes four CDC 160 words to equal one 1604 word, it was convenient 
to specify four quantities in the identifier. 

1. Current run number-indexes once a run 

2. Spare-may be sefit to anything 

3. Current block number in the run 

4. The number of blocks per run 
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For example, at the 5th run, using two blocks per run, and the "spare" 

set at 4321, the identifiers would appear in the 1604 as: 

0005 4321 0001 0002 
0005 4321 0002 0002 

Precise records must be kept of the identifiers when digitizing data, 
as data location on tape is by means of these identifiers. 

3. mm RECALL 

The data in digital form may be recalled by calling Subroutine DATA. 

S/R DATA in turn utilizes two subroutines, FINDIT, which locates the data on 
tape, and UNPACK, a machine language routine which unpacks the samples. 

The result is to create the original CDC 160 sample list, word forwor’d in 
the 1604. 

The upt>acked data is presented to the main program in integer format 
and appears in either KDATA (M, 1) or KDAm (M, 2) as specified by KLIST. 

Integer format output was chosen for the sake of generality as this Is 
the format of the raw unpacked data. 

Normal use is to change the data to floating point format, remove the 
mean, and divide all readingsi^ya factor of -409. 6. The letter factor converts . 
all readings to volts. 

When reading the tape on the CDC 1604, it is convenient to have the 
blocks arranged sequentially on tape to prevent tape rewinding and scanning. 
However, when digitizing data for cross -correlation, it is convenient to 
digitize all the right channel information, and then the left. 

r* 

Since it was desired to cross -correlate the left with the right channel. 
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some reordering of data on tape was necessary. Program SHUFFLE does this, 
reordering the data so that the output appears as right samples, left samples, 
right samples, etc. , in consecutive order. 

Once the data has been written on tape in the order that it will be re- 
called, the average time taken to recall one block, uppack it, remove the 
mean, and convert it to floating point form is less than 2.75 seconds. 

4. TESTING THE SYSTEM 

In order to test the capabilities of the digitizing system, a test tape 
was prepared with two blocks per run, and each block of 4000 samples on 
each of six known signals. 

Run • Signal Identifiers 



1 


8. Ov. peak-peak, 250 cps 
square wave 


0001 

0001 


0000 

0000 


0001 

0002 


0002 

0002 


2 


5. 0 volt, 250 cps 
Triangular wave 


0002 

0002 


0000 

0000 


0001 

0002 


0002 

0002 


3 


6 volt, 250 cps sine 
wave 


0003 

0003 


0000 

0000 


0001 

0002 


0002 

0002 


4 


-8 volts DC. 

(-3 below reference) 


0004 

0004 


0000 

0000 


0001 

0002 


0002 

0002 


5 


20 cps band-limited 
white noise, 7. 07 
volts rms. 


0005 

0005 


0000 

0000 


0001 

0002 


0002 

0002 


6 


White noise 7.07 volts 
rms. 


0006 

0005 


0000 

0000 


0001 

0002 


0002 

0002 



The data was analyzed by Program TEST and the resulting analysis is 



tabulated in Table IV. The data was recorded on Reserve Tape 18. 
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Table IV. Summary of Test Results. 



Run 


Signal 


Type 


|M?AN/EXP 

MEAN 


(J2/eXP^ 


fc/EXP 


1 


Square Wave 


250 cps 


0/-. 132 


16/16.365 


— 


2 


Triangular 


250 cps 


0/+. 086 


2. 0/2. 09 


— 


3 


Sine Wave 


250 cps 


0/+. 027 


4. 5/4. 56 


— 


4 


-8 Volts 


D.C. 


-3/-3.025 


'O. 0/4.74x10"® 


— 


5 


20 cps White 
Noise 


Bandlimited 


0/+. 147 


-/2. 10 


20/29. 5 
cps 


6 


White Noise 




0/-. 094 


-/2.2 


200/199 

cps 



Mean Recovery Time Per Block ^ 2.75 Secs. 
Digitizing Noise is 63 db Below a 3 Volt Signal 



5o ERROR. ANALYSIS 



In the analysis of the data it was desired to determine an estimate of 
the data second moment, and to determine a confidence interval on the 
.estimator. 

It has been shown QzJ that; 



s2 = 



for^ s^ may be conveniently calculated by: 



is an unbiased estimator 






s2 = 



N - 1 



- X £ Xj 



•ir I 






To obtain an estimate of the confidence we may place on s^, we assume 
that in the large sample case, as here, the distribution of s^ is normal in 
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which case; 




where k is the value of the standard normal r. v. which will have a proba- 
bility of being exceeded of ^ . 

The confidence interval here was chosen to be 99% or .99 which implies 



= 2.57 

ie. Pr N(0_jl) < 2.57j = 0.99 

Thus the confidence range C was constructed: 




C indicates the region about s^ in which one has a probabilily greater 
than .99 of finding the true^^. 

If we form C/s 2 we have 



"■^ 5 - X 100% 
s^ 




2 

N -1 



T 



X 100% 



In the data analysed here k = 2. 57, N = 3900 and hence C/s^ = 

<V 

5. 42%. 

This means that we are confident that the true^^ will lie within 5. 4% 
of the estimated<f2 = s^ with probability .99. 

Note that lim C _ ^ 
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In each of the runs it was desired to plot both the signal and its auto- 
correlation function. The units of the signal are con vertedlo volts by divid- 
ing all data samples by 409.6 or levels , 

5 volts 

Sine, square, and triangular waves were analysed to check the scal- 
ing accuracy. It was found that the accuracy of measurements exceeded 
that of the oscilloscope used to measure the input waves. 

The D. C. signal was analysed to observe the amount of digitizing 
noise and hum introduced by the process. 

The maximum change in the digitized level was about 4 numbers or 
^ 100% = 0.2% of the full scale. Over the period of 0.8 seconds 
the mean square power was 4.74 x 10”^ watts and the rms. voltage was 
2. 17 X 10“^ volts. (A. change of one number -2. 44 m volts) 

Thus for a 3. 025 volt signal, the noise and hum is 62. 8 db below the 
signal. 

The presence of a strong secondary peak in R(T* ) at O' = 1/60^^ sec. , 

indicates the presence of a strong 60 cps componentdn the noise. Assuming 
1 

that at T = ~ sec. the noise is uncorrelated, then R(l/60) = 2 x 10~° or 
60 

66. 6 db. with respect to a 3. 025 volt signal. Thus the digitizing noise is 
3.3 db stronger than the hum. If hum could be completely eliminated, the 
noise power would then be -63.3 db below a 3. 025 signal or -53. 7 db below 
1 watt. 

The "white" noise and the band-limited white noise samples were 
taken primarily to check the whiteness of the noise generator. Measurement 
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of the point of the RCT) curve indicates the generator output would be flat 
within 3 db out to 200 cps. 

A considerable discrepancy between the rms voltage as measured with 
a VTVM, and the calculated value of R(0) and s^, was observed. Some of 
the error is due to the motion of the needle, and some due to the fact that 
the averaging time of the meter was about twice that of the correlater. How- 
ever the most signifigant factor was a faulty measurement probe. 

6. RECOMMENDATIONS 

No attempt was made to "fix" the installation, and hence some prob- 
lems arose that would not have occurred in a permanent installation. 

Noise and hum reduction was only achieved by careful grounding and 
the use of BNC type coaxial signal leads. The use of coaxial cable is re- 
commended for all digitizing work, for the above reason. 

A clock pulse was used to synchronize left and right channel data. 
Unfortunately, the act of switching Ampex tape recorders on or off in the 
record mode causes transients to be generated on all tracks. This is a 
characteristic observed on both the Ampex CPIOO unit of E.E. department 
and the Ampex FR-100 unit of the Aeronautical Engineering Department. 

The transient pulses may be treated as clock pulses by the program, 
causing loss of channel time synchronization. 

The problem was partially overcome by removing the output of the 
pulse generator until the transients were past. This is satisfactory unless 
the first clock pulse is close to the transients. A standard interval of 20 
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seconds is recommended between the start of data and the first clock pulse, 
as a temporary restriction pending hardware modifications. 

A more permanent solution would be to imcorporate a time delay in all 
record channels of the Ampex CPIOO to allow transients to die out, when 
starting and stopping the unit, before energizing the amplifiers. This should 
prevent spurious clock pulses from being formed. 

Despite great care, occasional spurious pulses appear on the clock 
channel, and may cause false sampling to take place. The use of the short 
optimal length pseudo-random code as a clock pulse, together with a match- 
ed receiving filter would correct this problem. 

The system as implemented has a 2. 5 Kc upper frequency and 4000 
samples continous sampling capability. Since at best, samples may be taken 
every 190 jMsecs., of which only about 50 jAisecs. are chargeable to the A/D 
converter, the most likely place to begin the search for higher sampling rates 
is in the computer itself, rather than in a more expensive A/D converter. By 
deleting certain functions from program DIGITIZE, the sampling may be made 
faster but in no case is it possible to sample at intervals smaller than about 
100 jilsecs. 
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Table 1(a). Program Digitize. 



Cell 


Contents 


Code 


Explanations 


0000 


0101 


PTA 




1 


0603 


ADN03 




2 


7064 


JPI64 


Jump to INITIAL., 


3 


7500 


EXE.00 




4 


2410 


2410 


Set Enable 


5 


7500 


EXFOO 


1 


6 


1401 


1401 


Call A/D Channel 1 


7 


7600 


INA 


INPUT 


0010 


4176 


8117 6 


Store sample in- (007 6) 


11 


2076 


LDD76 




12 


-3464 3 


SBN65 


Enough Samples, yet? 


13 


6134 


NZF34 


If not, go to cell 0047 


14 


7500 


EXFOO 




15 


2400 


2400 


Clear Enable 


16 


2074 


LDD74 


Load Current Rug No. 


17 


4160 


STI60 


Store in cell 0133 


0020 


2067 


LDD67 


Load Spare ID. , 


21 


4161 


STI61 


Store in cell 0134 


22 


2075 


LDD75 


Load Current Block No. 


23 


4162 


STI62 


Store in cell 0135 


24 


2073 


LDD73 


Load No. of Block s/Run 


25 


4163 


STI63 


Store in cell 0136 


26 


7500 


EXFOO 




27 


2111 


2111 


Call 163 M.T.U. 


0030 


7303 


0UT03 


Output from 


31 


0000 


C 


computed L.W.A. 


32 


6102 


NZF02 


to i 


33 


0133 


0133 


0133 


34 


2246 


LDF46 


SetA-0137 


35 


4076 


STD76 


Reset running storage address 


36 


2075 


LDD75 


Enough Blocks yet? 


37 


3473 


SBD73 




0040 


6155 


NZF55 


If not, go to cell 0115 


41 


0401 


LDNOl 


if'-- 

Reset block no. ' 


42 


4075 


STD75 


43 


5455 


AOD55 


Update total black count 


44 


3454 


SBD54 


Check if capacity exceeded 


45 


6153 


NZF53 


If not, go to cell 0120 


46 


6061 


ZJF61 


If not, go to cell 0127 


47 


5476 


AOD76 


Begin intersample delay loop 
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Table 1(b). Program Digitize. 



Cell 


Contents 


Code 


Explanations 


0050 


2066 


LDD66 


Load intersample delay word 


51 


0701 


SBNOl 


If not zero, go back 1 


52 


6501 


NZBOl 


If zero, go to 0005 


53 


7056 


JPI56 


0760-M.T. U. block capacity 


54 


0760 


0760 


Total No. Blocks written 


55 


0000 


C 


address 


56 

57 


0005 

0003 


0005 

0003 


address 


0060 


0133 


0133 


Address of Current Run No. ID 


61 


0134 


0134 


Address of Spare ID of ID 


62 


0135 


0135 


Address of Cur. .Block No. ID 


63 


0136 


0135 


Address of Tot. .Blocks ID 


64 


0100 


0100 


Address of INITIAL 


65 


C 


C 


Address of last word of data 


66 


M 


M 


Intersample delay word 


67 


M 


M 


Spare ID anything 


0070 


M 


M 


Initial Run No. Set 1 


71 


M 


M 


No. runs desired 


72 


M 


M 


No. samples/block 76408 


73 


M 


M 


No. blocks/run . 


74 


C 


C 


Current run no. , 


75 


C 


C 


Current Block Number 


76 


C 


C 


Running Storage ^Address 


77 


C 


C 


JUMP CONTROL, 


0100 


4077 


STD77 


BEGIN INITIAL 


101 


2200 


LDCOO 


r 


102 


0137 


0137 


Set "A"-0137 


103 


4076 


STD76 


Initialize running storage address 


104 


3072 


ADD72 


Compute last address 


105 


4031 


STD31 


Store in 0031 


106 


0701 


SBNOl 




107 


4065 


STD65 


Store in 0065 


0110 


2070 


LDD70 




11 


4074 


STD74 


Initialize Run No. 


12 


0401 


LDNOl 




13 


4075 


STD75 


Initalize Block No. 


14 


7077 


JPI77 


END INITIAL 


15 


5475 


AOD75 


Update Current Block No. 


16 


5455 


AOD55 


Update Current Block Count 


17 


7057 


■ JPI57 


Go to 0003 



Table 1(c). Program Digitize 



Cell 


Contents 


Code 


Explanations 


0120 


2074 


LDD74 


Have enough runs 


21 


3471 


SBD71 


been done yet? . 


22 


6103 


NZF03 


If no, go to 0125 


23 


2074 


LDD74 


If yes, display last run no. 


24 


7701 


HLTOl 


Halt 01 


25 


5474 


AOD74 


Update current run no. 


26 


7057 


JPI57 


Go to 0003 


27 


7500 


EXFOO 


CallE.C.F. 


0130 


1111 


nil 


Write an ECF Mark 


31 


2055 


LDD55 


Load block count 


32 


7702 


HLT02 


Halt 02 


33 


C 


C 


RUN #ID 


34 


C 


C 


SPARE #ID 


35 


C 


C 


BLCCK #ID 


36 


C 


C 


BLCCKS/RUN ID. 


0137 


D 


D 


Data Storage 


> 


D 


D 




7776 


D 


D 


Data Storage 



Notes : 

1) C is ENTRY PROGRAM COMPUTED 
M is MANUAL ENTRY 

D is DATA STORAGE 

2) CDC 1604 Identifier will appear as follows: 



Word 1 


RUN # 


Spare 


Block # 


. Blocks/RUN Format 016 


Word 2 


Sample 1 


Sample 2 


Sample 3 


Sample 4 




etc. 
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Table II. Modifications to Digitize 



Mod 


Cell # 


From 


To 


Effect of Modification 


1 


0126 


7057 


7703 


Individual blocks clocked at end 



of a run. To restart run at 0003 



1 ^ 



2 


0014 

0015 
0117 
0126 


7500 

2400 

7057 

7057 


0300 

0300 

7056 

7703 


Initial block .clocked. Subsequent 
blocks taken as fast as possible. 
HLT in 0126. When run complete 
reset FFI manually 


3 


0014 


7500 


0300 


Initial block clocked. Subsequent 




0015 


7500 


0300 


blocks taken; immediately. All 




0040 


6155 


6113 


blocks have same ID. HLT when 




0126 


7057 


7703 


run complete in 0126. FFI manual- 
ly reset. 


4 


0013 


6134 


6506 


•V 

Eliminates variable sampling delay. 
Samples possible in this mode 
every 120 ^sec. 
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Table III. Program TEST 160 



Note: M means Manual Entry 



Cell 


Contents 


Code 


Explanations 


0000 


7500 


EXFOO 


Begin Bias -set 


1 


2410 


2410 


Set Enable 


2 


7500 


EXFOO 




3 


1401 


MtO/Vc/ 


Call A/D Channel 1 


4 


7600 


INA 


Input to "A" in Reg. 


5 


4070 


STD70 


Store "A" in 0070 


6 


7500 


EXFOO 




7 


2401 


2401 


Call D/A Channel 1 


0010 


7303 


OUT03 


Output Cell 0070 to 


11 


0072 


0072 


Cell 0071 INC. 


12 


6102 


NZF02 




13 


0070 


0070 




14 


0400 


LDNOO 


Load zero’s. 


15 


4071 


STD71 


Store in 0071 


16 


6414 


ZJB14 


Jump back to 0000 


17 


7700 


HLTOO 


HLT. End Bias-set 


0020 


7500 


EXFOO 


Begin Timer 


21 


2410 


2410 


Set Enable , 


22 


7500 


EXFOO 


Call A/D Ch. 1 . 


23 


1401 


1401 


V 


24 


7600 


INA 




25 


4143 


STM 3 


Store Sample in 0043 


26 


0300 


NOP 


No Operation-NOP 


27 


0300 


NOP 


These are Time 


0030 


0300 


NOP 


Dummies to 


31 


0300 


NOP 


Match Program DIGITIZE 


32 


0300 


NOP 


Timing 


33 


0300 


NOP 


... 


34 


0300 


NOP 




35 


0300 


NOP 




36 


2042 


LDD42 


Load Timing Control Word 


37 


0701 


SBNOl 


Subtract 1 


0040 


6501 


NXBOl 


If not zero, go back 1 


41 


7044 


JPI44 


If is zero, jump.to 0022 


42 


M 


M 


Timing Control Word 


43 


0043 


0043 




44 


0022 


0022 




45 


7701 


HLTOl 


End Timer 



53 























LJ 
































a 




























CO 




< 




























c 




X 




























CO 
































CO 




X 




























LU 




X 
























U_ 








X 
























o 




o 
































X 




20 
























LO 








C LNi 
















oc 


• 




O 




o 








o 














lu 






1 




II 






7T 


X 














-.0 


10 




to 








u 


w 


<< 














>" 


< 








Ll 




o 


< 
















5 






X 


o 


»— t 




X 


X 






















u 


1— 






cn 


< 


•x 
















LO 






e 


» 






cz 


— 




X 










o 


< 




X 


* — 1 


<T 


m 


o 


v: 


f\J 




tG 












O' 








1— 


\/ 


U 1 








1— 










►-H 




• 




o 


<I 




c< 


ry 


o 




X 


h“ 








>- 


lu 




(0 


UJ 


X 


o 


•— i 


L") 


o 














IL 


r'i 


i ) 


1:1 


V 




X 


CG 




f > 






G 










<r 


»■) 


o 


L.; 


>- 


ai 


U.I 


. — 


-J- 




>- 










h- 


t— 


X 


►-H 


‘.f 


CD 




o 


1 — 1 


— 




gg 


1— ^ 












TJ 




G_ 




o 




«N 


< 






X 








UJ 




\ 


LU 


X 


O 


x 


lU 




X 




n 










o 


(Ij 


c/) 


O 


X 


LU 


1 — 1 


X 


■— 


< 




UJ 


CO 








1—1 


X 


LiJ 






cr 


u. 


X 


<T 


G 




u 


to 












X 


n 


U) 


1 — I 








V 




X 










(— 




a 




GJ 


< 


O 


to < 






tu 


G 


e.'- 






»— 1 










ro 


1 — 


o 


cz 


X* 




01 


X 








o 


<r 


< 




o 






2T 


lU 


o 




LU 


< 


-j 






1 — 1 


uT., 


to 


PU 


h- 


o 


< 


1 — ( 








Ll 




LL 


«■> 




o 


o 








<r 


h- 


X 


X 






UJ 


V 


V 


O 






o 


Ll 


or 


<r 


X 


<r 




U' 


X 




o: 


u 




c 




-J 


cc 


o 


C 


I— 


Ll 


o, 


X 


X 


lU 






< 


H-- 


c 




<t 


a 






< 








X 


1> 


• 


o ; 


y 


c/> 






1— 




CY 


T— 1 


o 


CO 


X 






>— 1 


Z' 




X 


1 — 1 


— ■ 




u 


o 


UJ 






O 


X- 


10 


UJ 


G 


o 




0 


-J 






o 


vO 


m 


X 


UJ 


CO 




LlI 






7^ 


o 






K- 






1 — 1 




LlJ 


X 


CO 


Cl 


X o 


X 


rr 


X 


01 


« 


< 




vO 




X 


X 


L- 


X 


X 


X 


X 


< 


O 


<. 


to 


X 


o 




r-i 


i±j 


IX 


i — 






tCj 


<!' 


X" 


X 


u 






<c 


V 






X 






l— 


X 


X 


u 


►—I 


L_0 




1 — 






*> 




< 


\— 


UJ 


LU 


X 


< 








O 


<r 


to 


o 


•> 










X 




< 




u 


LO 


X 


01. 


X 


K-1 






r-i 




t/) 


> 


h- 


00 


X 


or. 


1 — 1 


t — 1 




Cl 


<■ 


X 


LG 




c 






OD 




1 — 1 




1—. 










c. 




X 


U.i 


c 








iO 




X 




r-H 


< 


to 


kD 


v: 


X 


X 


o 




O 1- 


Q 




1— 


c 


(X 




X 




X 




X 


<r 


t 


— 


— 1 


:z: 


UJ 




</) 


X 


< 


X 


< 


L-' 


1 — 1 


u 


K 


u 






CO 


ll) 


»— 1 


X 


1 — 1 




X 


1 — 


Cl 


<r 


X 


X 


h- 




<L 




1 — 1 




-J 


< 


X 


< 


u 






Cl 


X 


< 


X 


< 


H- 


^ , 




1-H 


X 






X 


vr 


X 


X 


X 


<1. 


X 




t~ 


< 


_J 


0\ 




fL 






< 






1.: 


X 


LJ 


X 




< 


c- 


a 


< 




X 






c 




X 


•—1 






IJ. 


ll! 


G 




•— < 


\— 


• 


to 


• 


• 




• 


''I 


X 


o 


UJ 


c 


Tl. 




lj_' 




< 


r-l 




f\J 


o^ 




<r 


Lu 




X 


X 




X 


X 




Z" 


n 














»— < 


X 


< 


X 






X 


1 


o 


V o 




o 


O 




X) 


LO 


X 








U) 




1— 


1— 




ai 




cr 


CO 




01 


c 


< 




X 




• — I 


X 


> 


to 


r?' 


< 




«.:! 


< 




<* 


X 


U‘ 




1 — 1 






X 


O 


S' 


O 
























< 


o. 


01 


Ll. 


:>2 
























X 


CO 


lC- 




























<r 


CCj 


Z) 




o 
























d 


G 


tr. 


c 


u 
























v: 


LG 





00 


































iU 


































G 


































< 


































uj 


































X 


































X 


































< 
































G 


































CO 


lu 




vO 








o 




















o 






r-H 








C\l 






















1 — 1 




o 








00 






















X 




•> 












X 
















O 


G 












o 




LG 




' 












vO 


o 




X 








X 




1 — ( 
















1 — i 


CO 




X 












X 
















\ 


X 












o 




NO 
















G 


G 




G 








o 




•> 
















CO 


LG 




U ' 












1 — i 
















O 






G 








LG 




— 


















X 




<x 








f— ^ 




<r 


















1— H 




LL 








II 




X 
















<1- 


G 




X 








U9 




< 
















O 


X 












< 




G 
















vG 


»— < 




< 








X 




















r-H 


X 




X 








X 




rv 






















<r 








:i/ 




X 
















X 


X 




G 








LA 




to 


X 














X 


1 — 








— 








1 — ( 


X 














o 






C: 




uD 




X 




X 


Ll* 
















G 




UJ 




< 








v: 


G 




— 










to 


LlI 








X 




X 






1 1 




X 






X 




X' 


00 




UJ 




X 




G 




X 






LG 






lG 




LG 


00 




< 




1 — ( 




►—1 




< 






t — 1 






1 — 1 






G 




Cl 












7-~ 






X 






X 




X 


U 






vT) 


X 


f\. 


LP 




G 


fO 




V 






V 




X 


UJ 




G 


f“H 


<r 




(X! 




•> 


fNJ 










•N 






o 






C 




oc 


00 




v: 


CO 










X 




UJ 






Ll 


«- 


r- 


•' 






LJ 






— 










CO 


DC 




7> 




X 


r-. 


X 




O 


X 




< 




X 


< 




UJ 


O 






c< 


2i 


r. 


X 




X 


X 




X 




< 


X 




r 


X 


. — . 


X 




Ij. 


oc 


*— . 




m 


t— . 




< 






< 




X 


01 


vO 






G 




01 






CO 




G 






G 






X 


r-^ 


►— < 




*— 1 




Cl 




— 


Cl 




X. 




•> 


VC 








o 








C\J 




r-H 






<r 






G 






L.; 


1 — 1 




o 


4T 


X 


00 


tA 


+ 


X 


iO 


•N 


VT 


c<> 








<_ 


X 


X 


X 


X 


«— i 






<r 


u 




t i 


•> 




II 








00 


e 


p'. 


1— ' 


G 






\ 


<c 


c 


II 


e 


1 




cr 








LO 


P'. 




X 


tr 


(X 


X 


CL 




7>" 


G 




Z 


c<^ 


UJ 


X 


— 






-- 


•— t 


< 


CO 


< 


X 


il 




OC 


X 




CO 


G 




X 


X 


X 


X 


Ll 


X 




1_ 






O 




< X 




X 


CO 


< 


< 


<L 


<!' 




X 


X 


11 




UQ 


CM 


X 


X 


c\i 


X 


t-H 


X 


li. 


O' 


z 


lo 


X 


»— 1 


X 


X 


X 


< 


00‘ 


X 




00 


X 


X 


X 


CO 


01 


01 


if 


X 


— ' 


1 — 1 


< 


X 


X 




1 — 1 


II 




1 — ( 


X 


X 


o 


c 


o 


C 


< 


u. 


CO 


z~ 


<c 


IL 


O 


CO 




O 


X 


c 


< 


X 


X 


X 


X 


UJ 




X 


G 


u 


v: G 


X 


G 


G 


X 


UJ 



o 

«cT 

Li. 





lO 


00 


c 




CC'i 


vO 


X 


Ch 


CM 


CM 


CM 


c<^ 


CM 


(MJ 


CM 


CM 


PJ 


00 


CO 


CO 


00 


vX) 


CO 


00 


00 


OC 



LJ U 



U O U L./ 



U U U U U 



54 

















































Llj 








































O- 








































G- 








































O 








































h- 








































c^} 


















f-H 






















f~ 








































O 


















CO 








































o 






















to 


















( — 






















<j; 








































"T; 


















O 
























vO 


vO 






































r-J 


♦—1 














.r, 






















ro 


o 


c 














C' 






















rG 


















II 








































tS* 






















h- 


CC 
















<vL 






















D 


llJ 
















X 






















CO 


h 


X 














Ll. 
























<* 


r.P 














h-~, 




lf> 


















o 


i.' * 








— 
















1 — i 














u.f 


















\ — 




OD 




Ovj 














i — 




' 






J — 1 
















CO 










— 




u 


x 


cc. 






X 








1,1 




C 














r-< 




_Li' 


<0 


< 


— 


CL 


< 








o 




i — 




O 










C: 








IC 


O 


X 










1— « 








X 








c 


Cv 




Llj 


o 


IT 


r—i 




•— 












r 












c' 


r~i 




n 


LlJ 




o 


“0 










a; 




C 




C' 








-J 








o: 




€\ 




LO 








r-l 


tr 






e* 








li 


V 




to 


Cr 


X 




o 


C 




c 




a 


cz 


L 




















_) 


< 




X 


_J 




uO 




X 


0*1 






1/7 










C 




■*r-' 


u 


o 


LU 




'Q 




r> 








a. 












X 












i — 


r— ! 


1— t 




O’*' 




t— * 


0- 


> 




cj 








< 


CO 






o 


o 


< 


r-i 




o 


rH 




or 




X 


















c 




Ll-' 




CO 




r— ' 


CO 




d. 


CO 


X 




<1 


X 






<Tv 


»> 




Ol 


CC 


I — 


O 




f— H 


CO 






t-n 










il 






f"" 






or 


o 


< 


>J 




— 


•N 


Lf^ 




X 


i: 


c 




0 


to 








CX) 




U.‘ 


or 






LO 


xc 


CC: 


o 




X 


C' 


i~ 




<*. 


-iXi' 






L’ 








iC 


o 


o 




KJ 


O 


03 




u 


CC- 






X 


-J 






n 


c 




X 


LiJ 




X 


iO 


o 


00 




— 


C' 




X' 




Ll 


Ll 






•— 


( 


V .- 


X 








uO 








X 


r-^ 


— 


r— i 


O' 


»— < 


t — < 








c 


1 ,' 


1— , 




i-i-’ 


ll I 


< 


•0 


X 


— 


r— ; 


•' 




00 


r— ( 










h- 




c 




h- 


X 


■O 




K— < 


c 


1— 1 


or 


ii' 






00 




X 










_J 


c 


« — ^ 


< 


00 


CO 


— - 


a; 






f — 


— 


2' 












9 


<■ 


C’' 


o_ 


cr 


X 


< 












vi: 


V 


C.''’ 


C' 


^— 


X 








h- 






<j. 






r-H 


1 — i 


vO 


L,' 


r-H 


|.— 


'v. . 




1- 










f— . 


< 


»> 




Q. 


t— . 


b 






C-. 


Q 


CX) 




C “ 


t/: 




Ll> 


LlJ 






Li. 


r 


<r 










-r 


r-i 


CO 











ic 


O 


G: 










V. 


• — 


c. 


(T 


r- 


k. 








it :■ 


«> 




0, 


<r 


ts 


K-. 


V-l 






i 1 ' 




<■ 


‘.’!- 


-.1.. 


~j~ 


•U_ 


-j' 






1—. 


CL 


0- 


M- 








• > 






■ 


~ 


r_' 


<!■ 


•V) 




r-: 


\ 


TG 


r-; 


1 




( '- 


! 


ti.. 


X- 


CO 


CO 






»— < 


O 




'JT- 


CO 




:\i 


X 




r' 


■ — 




X 


X 




1 — 1 


r— i 


f\' 


Ll- 




\~ 








— 




— 


< 




• — 




b 


r— t 


irtt 






a) 


CO 


X 




3 


fj) 




} — 


! — 


1'- 


; — 


'r 


cr 


1— , 


lT: 




<X' 


Ll 


X 


o 






7C 




c 




o 


<; 








U 


LiJ 


G’ 


Cl> 


o 


X 


C' 


<10 


LT 


X 


X 


i—i 




Q' 


lL 




X 




* 1 


* *_ 


r-1 


i-i- 


Zj 




1 — 


X 


►— 




»— 1 


Zc 


2T 


X 




cr 




I 


C< 


or 


0; 


•V 


X 


ii.. 




— 




,-H 


-- 






1— 1 


V—. 




0 






c 


o 


o 






<. 




L.- 


Ll 


c 


CO 


Ll. 




Ll.. 


CL 


or 


0 


2 : 


U. 


Cl 




U. 


Ll_ 




Li- 




'0 








Cl 




or 


£V 


X 


Q- 


0 


Lii 








vT 










iC\ 


V O 


f'- 


CIJ 


C 


C 


r— ‘ 


CM 


r' 


o^ 


r-l 










f— ’ 






Xi 




o 


c. 




O 


r'\ 


r-i 


1 — < 




r-i 


f— i 


(M 










CC 


cr 




O 




CO 


CO 


CO 


00 


CC 


CO 


CC 


CO 


a) 


CO 


00 





55 











































u 






















(/) 




-;:r. 














< 








o 




■fYl 






— 






<T 




















h- 




1— 




< 


■«— H 










OT, 






0“' 




•J) 


<r 


O 


1 
















l-H 






o 




X 


lO 




o 




b 






-J 




1 


V 


fv 


< 


o 




o 






1— 


O 






V 




t/) 




CV 




< — 1 


00 


^- 




(V 


*N 






H- 




”> 


< 






cy 


u 




< 






C\J 


U1 


X 


II 






< 


o 


— > 


o 




— 






>~4 


i*n 


1— 






1— 


<r 


O' 


u 


f- 


<r 




c 


— 1 






re 




<r 


<r 


Kf 




X 


h- 




H* 




n 


Q 






o 




Z! 




Li I 


< 








U : 


"Z~ 


_j 






r_ 


1— 1 


O 




c 




rj 


X 






-j 


r>^ 


z: 


►-H 




G 




L/ 






<< 








o 


►-H 




CTi 


— ! 


U- 










< 


r— 1 


u 






V 






O 


{ U 




Li. 


“S 




+ 






o 


u 


> 






o 




O 








Lu 




tJ_i 


o 


(Ti 


►-H 


lO 




ty- 




>/ 




X 


o 


X 


LZ 


— 1 






re 


oo 


re 


oO 


u 


KJ 


X 




<r 


u 




Z 


X 


o 


ry' 


o 


fV 


o 


a 


< 


q: 








►—1 


< 




G 


< 


b 


-J 


Cl 


*T" 


o 


“0 


a 






*>. 




<T 




< 


m 


re 




o 






u_ 


c. 








Li* 




►-H 


IT' 


X 


< 


U1 


V 


o 


L 




LlI 


Ll' 




LiJ 






< 






u 




— ! 


Ll. 


O 


re 


<rr' 


•re 


1/) 


ce 




00 




o 


L/> 


-J 


re 




o 


< 


o 


o 


CC’ 


s 


1— 


rH 


— 1 


or 




C 


C’ 


^— 


> 


h- 


CL' 


3 




U.J 


r*) 


CO 


o 


Li_ 


K 


o 


or^ 


O 


(O 


<; 


lO + 


-h 


o 


+ 


»— < 


+ < 


+ 


Kjy 


+ 


■r 


<■ 


+ 















< 






















V 


1— 












— 










u 


< 




LO 




h- 




f— 


\— 








O 


O 




LlJ 




oO 






U1 








.J 






re 


o 






►-H 










(T! 


00 






X 


_J 


e 


J 




00 






►-H 


►-H 




< 


LlJ 


\e 








ce 














y- 


<rs 


o 






o 


L- 




ZT 


< 




o 


X 


<J- 


y- 






<r 


L/) 




►-H 






\— 


U) 


+ 


z 














<r 




z 




H-H 


,-H 


< 






►-H 




G 


o 






z 


— 




1“ 


R- 


t—i 






oe 


V 






‘G> 




o 


< 


o* 


X 


L- 




o 






o 


►-H 




ce 


O 


b 


LiJ 


I— H 






z: 




ce 


oO 


L- 


o 


V 


X 




X 


c\i 




>— 1 




o 


9\ 


< 








h- 


l:J 




c 






5i 


L— 


o 




or 


rv 


LJ 




X 


f— 


o 






lO 




o 


o 


b 


1 jJ 


S' 


iJJ 


<: 


ce 


■0 


X 


X 




re 


Li- 


Ll_ 


ce 




.o o 


o 


— 


1— 


X 


z 


o-^ 






ce 




re 






V 


<r 








o 


C) 


o 


IL’ 




1 — 




KJ 




1— 




y— 


1“ 


1— 


u 


nr 




iO 


h- 


O 


1— 




<1 1 


u. 


L- 


h- 




O 


LiJ 


ce 




-J 


Lu 


o 


re 




h-H 


h-H 


1— 


1— 


> 


►—1 


ce 


rn 


►-H 


h-H 


o 


X 


CO 


o 


LU 




< 


Li. 




H-H 


X 


cc 


y- 


oO 


+ 


+ 


00 


T 


oO 


•i- 


LU 


+ 


oO 


+ 


co 


+ 



vf ^ v: 



u 


‘u 




U 






g; 


u 








+ 
















< 


!~ 


<r 






ce 


O' 








f— 


i — 












a. 


a 


►-H 


G 






X 


c: 








r-, 
















z: 


X 


Z 




rH 


r— 1 <■ >:<: 


<■ 






< — 1 


X 


X 






c\J 


sD 






— ) D' 


LU 


X 


G 


X 


(Ni 1 <r 


cr 


-H X 


1 — I ”0 ' — i ”0 ' — 1 


X 


Pm La' 


LU C\. 






rH 


CO 








rH 


rH 


rH 








rH 






rH t— t 


C\J rH 




1 — 1 






C\i 



y~ 

z 

111 


X 

ce 

X 


X 


X 


X 


< X 


oO 


<n' 


. 






< 


. 


«< 


X 




X 


<r 


X 




X 


X 


^ r-H Z 


X 


oO 


X 


o 


z 


X 




Q < 


ce 




< 


s 


< 


X 


<■ 


-y 


< 


“V 


< 


z 


< 


z 


r— 


,—4 


z z o 


fV 


re 


1 — 




Li-' 


00 


oi 


X 


X oO 


<r 


►-H 


CTj 


X 


L/) 


X 


00 




oO 


►-H 


oi 




oO 




i/v 


on 


u.. UJ X 


X 


o 


CO 



V/ 

u 

<c 

a. 



Of*) 

O' 

c. 

< 



tNJ C<“ 



56 



QRS 36 RIGHT JUST ,SIGN EXTEND 

STQ 2 +STORE IN KDATA ( I +3 , KL I ST ) 

LRS 1? SHIRT 2ND WORD INTO A REG 

ORS 35 + RIGHT JUST, SIGN EXTEND 

STQ 2 ** STOR- IN KDATA ( I + ? , KL I ST ) 



o 

CO 

H- 



Z) 



• 






a 








V' 


— 






o 


r-H 






U 








c 




C\,' 




*-< 










X 






a 










LJ 






ZL' 


c 




X 


1 — 


O 


LiJ 






Z’’ 




<c 


<.i; 






















5 — 


L!J 


r 




-0 


a 






z 


ar 


.“V 






u 


Li 1 




o 




6 




z: 


ce 


(z: 


LI.' 




f— 










c 




Oi. 


”0 








f— i 


h“ 


O 




c 


f— 


4- 


V/ 


4- 


O'. 


f— 




a 


U'. 


►— i 


C.0 


"0 


u; 




*n 


S 


1 — » 


1! 


1 — t 


n 


C*' 


Li.l 




3 


+ 




+ 


"0 


4- 




4- 


O) 



kr 

Or 

C' 

■,% < :,;c 

'I" 



C\j (Nj r-i r-i rv 



< ^ V "Z — c 

h- :<: oo _j .v: z: _j z* 

^ tl, ^/) Uj 



<r 



C.0 

ri h- 

<•■ X 

re l: 



57 









lu 










































u 






































CO; 






c 
























JZ 


























h- 












Q_ 




• 


Cl 






i Ll 


















UJ 






































c 


u'- 


o 


L.' 






f — 


















o 


























h- 












c 


LX 


.-IJ 


C 






<l 


















f— . 




































c; 


lT, 


LJ 


"V 


X 


Li 


) 


CiC 


















lu 












1— . 














# 










a 


X 








U 






















xi 








X 




— 
























c.^ 




c. 


cO 


u 




- 


CO 


>- 
















o 








< 




X 














fXf 












X 


lO 


h- 


ex 


c 






(X 
















LJ 








2 “ 




h- 














»— 










c 


!— 


X 






"V 


- 


►^-1 


lU 
























n 




X 
























t-ri 


d 




o 


Ll 








> 
















LlI 








cO' 




Lli 














Ovj 










CM 


o 




> 


u 


L.f 


i — 


a 


ll I 
















(X 








a 




,:Z 














h- 












h- 


LLI 


<1 


f— t 


f— 


U 






• 














< 








o 




t— . 




n 




















b : 


X 


X 




O 


r-^ 


(J 


«rr 


Ci 


U- 






























o 










Lf ' 












< 




LI 


O" 


X 


> 


C.0 


U! 
















o 














Li • 


o 










C'^j 










< 


c^ 


u 












i 


•'1 














CO 








*>■* 








-J 




























Ll 




1 -> 


CJ 


o. 


l- 






















X 




>— 


t— 4 






• 


















CL 


a 


LJ 


1- 


li' 


L- 


NZ 




c/) 














<r 








Z' 




1— I 


h- 


n 


< 


o 








o 










t 


i 


o 


1 — t 


Cl 




in 


< 


X 


















• 




Cl 




f— 




tO 


h— 


o 






X 


<o 








CO 


a 


ex 


> 


X 


— > 






1..0 


O 














CO 




X 




CO 




X 


X 




X 


c> 






f— 1 


1— i 














c^ 




a 








u 














1 




r-C 




O 




iiJ 


J 


^-0 


■-r 


c 
















n. 


iT\ 


vO 


1 


I 


( 


-H 


L^' 


lO 




• 












cO 




(U 




tO' 




n 


-T 




X 


* — 1 






r-H 










1 












— 


O 




UJ 


L- 
















<r 




f— i 




t— I 


U 


< 


II 








h- 










n 


c\. 


X 


X 


Pvl 


X 




_J 


« 




















X 












' ; 


X 


cj 
















;> 


C' 


c 


C. ’ 


c 


c 


'..I 


aj 


li. 
























f— 




X 






, r 


— 1 






Ij J 


CO 








c 


c 


c 


c 


e 


c 


*0 


\ 


</.) 


I— 












vO 










Z"’ 




o 




cO 


h- 








n 


« 








• 


c 


c 


c 


c 


L 


X 








►-H 
























CO 


U'’ 


1— . 


X 










— 








cr 












X 


Ll! 


c 


<r 












O 






• 




cr 




1 — 1 


1— 




< 


iC\ 


















X 


X 


X 


X 


Cvi 


. 1, 


-J 






X 












X 




X 




a. 






a 


rH 


' r 


II 








c. 








X 


c 


c 


c 


( . 


C„ 


: .' 


a 




Ll 


o 










h“ 


< 




r-H 








u- 


t— 1 


rH 


X 


m 






(M 


1 — < 








c 


c 


c 


c 


c 


c 


1. :' 


“ 


Ll ! 


C* 


Lu 












u; 




UJ 










X 






















i 


c 


c 


o 


c 


c 




< 


h“ 














fa 


Xl 












t— . 


cO 


• 


• 0 








c 










r 












O 


CO' 




X 


O 










.'Z 






>: 




r— * 




a 


O 


z 




LU 






c ■ 


< , 








X 


c 


c 


e 


c 


c 


\ 




X 


p 


LlJ 










*— » 






1 — 1 




C" 




a 


2^' 


u 


rH 


\ 






< ' 










( 


o 


c 


o 


0 


C- 


< 


o 


>: 




CO 




















C.' 






c 


< 










<} 


I.'. 








c 


c 


c 


L 


c 


c 




c 




C' 


I" 










1— 










r-l 






X 


ft 


‘O 


(X 






*— 


•—■ 








c 


C' 


c 


C 


c 


c 


f — 


c 


C 


1 — i 












<0 






• 










c 




h- 


\ 






'0 


c/) 








C‘ 












a 


'J- 


Lu 


\ 


U' 
















X 




f— 






o o 


lu 


CO 






1- 


• • 








c. 


X 


c<^> 


^0 


lO 


vO 


Uj 




h- 


T~i 


• J 










<c 






1 — 1 




X 






rH 


Ll! 


t— • 


I— i 






•'•'C 










c 


C:* 


c 


c 


<w 


f. 


> 


X 


1— . 














1— 






LU 




t — i 








h- 


X 


\ 






o 










C' 


c 


r.. 




o 


c. 


X 


00 






1— 










a 










cx 




X 


H- 


ZL 


CO 


o') 






X 


c 




— 




1—1 


c 


c 


c 


c 


c 


(0 


1 — • 


»— 1 


lC 


CJ) 










o 


a 




X 




a 




K 






CZ 


\ 








r-( 




< 




c. 












O 




_J 




Z.^ 












a 




r-C 




rA 




X 


►-< 


''y 


X 


o 












1— 




o 


C'. 


c. 


o 


c 


e 




UJ 


1 




X* 










X 


►— » 












LlJ 


iX 


(X 


r 


<}- 






CX' 


O' 








c 


X 


X 


X 


X 


X’ 


o 


a 


Cs 


X 












cO 


u. 




<r 




X 




O 


a 




X 


\ 






«• 


l- 




r' 








<-■ 




<< 




1 — 




X 


c. 














1 




• 












C' 


a 


o 






o 


ll 




X 




O' 














1— 


<r 




rH 






vO 






h- 




X 






CO 


C'. 




X 


■H 


\ 






cr. 










X 


X 


Tv, 


r-vj 


X 


X 


CO 




X.' 






— 




( — 1 






zr 




rH 




•1 


O' 


— 




<J. 




r-H 






r 






•-"T 




< 


r 


c 




r 


o 




U' 




“r; 




<r 




o 






L'J 




UJ 




LO 


c: 




CO 




II 


\ 






<r 


U: 




l“* 


c/) 






c 


r 


c 






> 


Ll 




h- 










X 


C 




r- 




a 


~r- 




h- 


fZ 


X 
















h- 


CvJ 


C' 


c 




c 


V 


V.O 


Cr: 


-L. 


O'”' 


1—1 






X 




O 






X 




Z) 


a: 


r-H 


L^ 


< 


<< 










«r 




c:'. 


ll- 


c 












• 


u*. 


h“ 


c. 




>c 




CO 




oc 






X 




z 


"T' 


h- 


4— t 


LV 


h“ 


;q 






1— 






V / 






1— i 




r-1 


r— i 


r— > 


CO 


cO 




X 


,o 


CO 








X 




— 


«' 




iZ' 


b 


X 


X' 


a 


X 








<r 


C' 




— 


X 




C' 


C' 


C' 


c 


c 


o 


Ll. 




o 








vO 




LU 






vO 






4ir' 


Ll 1 


CO 




< 


X 




i— 


r' 




< 






r— t 


V. 


V. 




r 


r 


<J* 


cx 


in 


u 


Ll 


m 




rH 










.— 








CZ 


O 


X 


IZ 


o 




CO 


X 


r— 


h- 


b- 


>/ 


{. 


c 


c 


v_. 


c 


c; 


1 






IL' 


a 






r'> 


1 — < 


f-H 


c 


< 


c. 




X 


C' 




X 


X 




CO 




LlJ 








U 




c 
















X 


Cj' 


c 


r. 


r— i 


*' 






X 


2 ■ 




CO 


c> 


r-t 


C'. 


C' 


< 


>o 






r~ 


X 




C 


X 




c 


c: 


C ' 


cc 


c 


c. 




6 


> 




p 


X 


X 


X 


X 


r-- 


oc 


X 


X 


< 




II 


c\j 


CO 






h~ 






c5 


<' 




LlJ 


‘•r 


L. 


c 


c 


o 


CO 


c. 


cO 




a: 


Cl' 




c<*^ 


r— i 


.o 


'.0 


X 


vO 


•— 




vO 


c: 




O 


cz 






h- 








h“ 






bJ 


e 


o 


c 


o 


c 


Cl 




CO) 




LJ 






--- 










uO 


— 


— 


cZ 




o 


o 


b 


< 


IJ i 




< 


c/' 


0' 




< 


or 


r 


c 


c 


c 


c 


c 


'0 




h- 


1—. 


< 


h- 


h“ 






. — 


h- 




1 — ’• 


h- 


i—i 




r-H 


i-H 


X 


X 


cy. 


« 


£T 


or 




c 


> 




C 












7j 




a 


X 


CJ 


<L 


' i 


X 


<c 


< 


<. 




•tc 


< 




o 








'll* 


c< 


Uf o 


Uj 








LJ 


r-l 


X 




<r 


lO 


vO 


1 — 


< 


Ll. 




H-H 




;s^ 




*■> 


u. 


X 




X 


X 


o 


1— 1 


O 


C) 


< 


X 


< 


•> 


C) 




X 






Q 


r 


e 


c 


c- 


c 


c 


LJ 


h- 


l; 


<L 


O 


ro 


rr 


o 


a. 


V 


a 




a 


O'. 


< 




‘Cl' 


<x 


1— 


<x 


'.'O 


<l 


rt" 




f-“ 


■'C 


6 


>- 


c 


O' 


c 


c 


c 


c 


<i. 


< 


X 




c 


C') 


O 


V.Z 


o 


O 


o 




c 


C 


UJ 


O 


LX 


u- 


< 


Ll 






Cl 


c- 


— 


u 


IL. 


f— 


c 


c 


c 


c- 


c 


c 


<i_ 




u 


a- 




Lu 


Lu 


lu 


IL 


b_ 


IL 


LU 


X 


Lu 


oz 


o> 


a 


a 


o 


X 


«\ 


— 1 






1 — 










































rH 


















a 








































i-H 


X 


rn 


<r 


lH 




vO 


u- 






0 ; 








o 






































O 




o 


' z 




o 






CO 






f-H 








o 


h- 




































o 


c 


C 


JJ 


o 


r2> 






C' 














u 

( 


u 

1 














c. ■ 


X’ 


LJ 


u 


L 


Ly 


U 


L.' 


LJ 


U 




rU 


1 — ' 


1 — i 


r— 1 


r-t 






r— 1 


1 — ( 










LJ 





58 



Ll 

o 











































vC 


















10 


• 




































»— t 














• 


















/ J 


r-H 




































1 






v/j 








O' 


















0 


i 




































<'.r 






• rT 




























■< 




































J 






J] 




















no 






'0 


<t 










• 


























L-J 






Ll 








1 












cr 






0 ; 








rv 


























f-H 






































X 




































cr 






•x 








'Xs. 












:q 








•— 






h- 




Cv 






















X 










f— 
























K-x 


rr 


X 










































'./) 








r- * 












1 




r.G 


fx’ 


— 






<7 




V/': 








































• X 












— - 




.0 


6 


CNj 






Jj . 






























1 — 






,J 








') 












i\j 




lT5 


OJ 


'•7 






Y 




. A 














• 








h- 




' 






\/ 




















•x 




X 


i 


“0 






C 






























' r 


























"n 






X 


•y) 






• 




o 














■ J 








J 










X 








1— 
















iJ.. 




r 






r—l 




c 






















X 




L‘_ 






'X 








. : 












rri; 






■ L 


nj 










i! 














1 


o 






(i 




t— 1 






TY 








(Y 












h- 




>1 




Ql 






X 




,—A 














— 








— 










r« 








X 




Xi 








-f 




— ) 




<x 






nj 




(.n 














X 








(\! 




i.C 






t— 






--: 


ii 
















X' 


'.0 


n 






y' 




•-• / 














■ — 








•' 




r- 






X 






«X 














X 




'Ii 


h— 


X 






if) 


















— 1 


r- , 






/ 




_j 




. ^ 


I'J 






“) 


r-l 


nj 














CO 


-I 


1 


























*\ 










Cj 




■ ? 


:.::i 








• ' 


•' 


pv 








:i 




<y 




(\i 






1 






























O- 












o 






r. 












■~r' 


> 


<n 






f\> 




• 


r-! 




























-> 








1— 


— 




—7 








^xj 




X 




'0 






if 




f 


It 












• X. 


< 






l~ 




O 






< 






* L 


<7. 


•Y 


— 








•' 






D 




f\i 




lO 




ti 


! — 












i — 








<i; 




h- 






»- 






C' 


i- 


H 


•Y 




X7 








I 


1 — 


II 






11 






U' 












' c 
















•V 


.--t' 








o* 


<r 


K- 




<T 












nj 


lO 




nj 








— 












n 






>r 




'.I- 












(I 




r:.- 


•rl 




!• 




< 




wl' 




LO 


— 


--- 


Ll 




IT' 


— i 


;'■ 










X 


X 




'y; 






h ' 




• ,o 










X 




c. 




X 




1 — 




L: 1 


U J 


<0 


lL 




P? 




cr 














!l 


















-J 






l“^ 




h' 


X 




v.^ 




<r 




l' 


1 “ 




h- 


— 


0 








_J 












Y 










(<' 




') 


_i 






• 




IX 


1 - 




n, 




0 




< 


''C 


(4-> 


CO 


(■'' 


1 . 




C' 










(■' 




r-. 






•' 


T--. 




L; 






<L 












fv 








> 






• 




c 


• 










'■ 










r 


«— 




f 


o 










OJ 






w 






i.r 












C) 


3> 


c 


ir 


r*— J 


L-"'- 






r 








•' 






( n 




ii 


o' 
















X 




K/^ 






LO 








0 ) 


X 


• 




i 






• 


— 


— 












( tj 










> 




:.) 








! — 


• 


11 


II 












*X 




1 " 1 


II 


X' 






(, • 


1 — ( 


u 




i 


r. 






It 










OJ 












- 1 ' 


( , 


(X 


n. 




11 




r~i 




r-l 




1 


nj 


< 






ii 


I 






— 1 


•' 




1 — 


r- : 






r 
















.'-j 


X .> 


' T 


Oj 








cr 




0 




X 


*rt. 




lY 




r— » 


>■ ' 


j 






(\: 










**1 


- — 




O 












X 


vX 




0 ' 




Cv> 




0 : 




..J 


n 






X 


0 






i' 


\/ 




o 


I'Z 






<.f 




f~> 


1 — -1 




-%• 






It 








i 


I?- 






70 




'■n 




c/) 




>■ 


o3 




1 — 




V 


Y. 






h- 


-- 




>c 












• t 




r-t 


1 — 






tA 


— 




LA 




<r 




X 




X 


r-r 


X 


r~. 


X 


0 




c7 




X 










iO 


Vr: 




O 


X 








iT' 










X 


if) 










1 










m 


c- 


< 






II 


' : 




. ') 












o 






n 




^•r 










— - 




*n., 




'y 




— 




Ll 


0 


X 




• 


U. 






K- 




r ■> 


' -"Z 






i—i 












1 • 












rv 


nj 




V 








r-l 




0 


111 


— 


»o 


fXi 


-f 








c> 






'•X 






r-l 






•c 




Z:» 




X 


,/ 






*■ 




1 — ( 






i ') 




rx 






>’ 


I — 1 




>— 


c: 






1 : ' 


1 — 


• 




1 — 










< r 


!— 




• Y 








r 








• 










n 




'^^ 1 - 


C) 


•T 




!i_ 


• 




; ■•'' 


o 


Z’" 


i-H 




;*•' 
















” 




!1 


•O 


tn 




X— 






r- . 




X 








;5 


'T 


— \ 


X 


H- 


r—l 


i“H 


i 


♦— 


!.i 


t." 


h~ 














c 




li ■ 








•X 








— • 


#X 




<T 




- 'X 




lO 


IlJ 


<0 


r—l 


ir 




a- 


in 




c 






C' 




1 — 


■Ct 






X 




Cy' 




V 






X 


t — 


f— 












s— 






iV. 


:Ic 


0 *} 


C3 


X 


rH 




i-f- 


*— 1 




ui 


y~^ 




• ; 


.-- 














Z'L 




tn 




• ,' 




h- 










0. 








1 — ( 


(y'» 


if) 


r-i 




’ '. 






: A 


CO 


I 












it 








< 




•X 


*>' 






..J, 


<r 




X 


X 


r~'^ 




r-l 




{Y 






lO 


U- 


“ ' 


t-i 








— 




• / 






> 






< 






(\) 




*• 


X/ 






1 — 




r-i 


<■' 


>c 




1 


< 


<l 


1 '. 


* 


m 


L*" 


“ 


C 






• 






ti 


X 




X 








' ■ 


X 


I- 


lf> 


rH 


II 


ll 


>L 


< 




'T 




ii 




, •’ 


Ui 


cn 


I— 


r- 


— 


0 . 




c. 


1 — 




•t 








H 




C-- 






' 














-- 


+ 


l’•“' 






•X 








> 


X 


r-/ 


m 


1 


a. 


r— 


1 — ! 


< 


<;■ 


,z 


«~ 


r* 


r-\ 


f— 






•N 






C'' 




lU 




il 


Osi 


(X 


I — 1 


X 




LT; 


f — 1 


. — ( 




X 




1 




• 


r— i 




il 




o 


_J 


•Z-. 


\- 


i. 


9' 


„ • 


Ll 




V/ 


L. 






C— ( 


o 


CD 






«X 




-f* 


u . 




1 ! 


«x 


1 ... 


1 -H 


Y: 


1 — t 


‘.0 


f\; 


fj' 


Y' 




1 — 




1 J. 


r-n 


n' 


— f 




“"'l 


•. 


(> 


\^/ 




r ) 




u" 


I—I 






■~) 


“) 


b 


r— 


ZCj 


II 


”:) 




■3 




0 


10 


11 


II 


1/; 


Lu 






_l 


V 










lY 


~2' 


ti 




XO' 


Y" 


Ll 


Xt 


il 


-J 




— 




uc 


u 


<L. 






— 




Ii 




uO 


r—l 


iV 


II 


0 


r-; 


*— ♦ 


_J 




r-- 




f\ 


-X 


ii 


t— . 


r—t 


<: 


I— . 




O 


II 




Ll. 


o'^ 


<r 


*x 


II 


u 


I—I 


r-H 


in 


<; 


I— . 










0 


r—l 


LO 


r-H 


O' 




il 


. i: 




t—{ 


1 


f— i 


1— 


o- 


h— 


l~ 


<L 




'XJ 




V 


f— H 


1— 


i— 


I — 1 


ii 


1 — 


fV 


r—l 


t— 


h- 


nvi 


< 


Ii 




h- 


U 


C.' 




CO 


VJ 






«— 




tr 




X 


<0 


<1 


7'' 


h- 


u 


uO 


0.1 






a 


<' 


n 


r--! 


X" 


<1. 








'X 


r— 


r-l 


0 


C.J 


70 




(7; 










i.L 


r , 


COj 


IT' 


r 

\ 


lT 


C- 


c 


<i 


< 




■ Xj 


lL 


L..* 


O 


c 




i.' 


C 


c 




CO 


c 


'I' 


< 


1/) 




< 


0 


r— , 










a. 






X 


V/. 




> 


X 




r:; 


-• 




1 — 1 




o 


X 


X 






V.; 


X 


0 


X 




lO 


0 


o-j 




IL. 


V.-* 














r-- 








CY. 






-1 








C\J 




<r 










o' 








lC\ 






m 


















S^i 


'. r> 


. f' 






f— 1 






r-^ 








in- 




tn 










i-H 








rH 






in 




















59 






c/) 

1 — \ 

li_ tXi 

ii- 

O 

O 

^ U 

r—i 

v/) fVl 

to L'^' 





vO 




















rH 














c\i 






— . 




rx 


















Xj 


•i: 


(X. 












3 








rfj 


<* 
















r-H 




X 


a’i 












•-'1 




1 






X 












h- 




3 


.-r 




















H 


f — 
















cS 




— 


< 
















X 






0 


1—1 


W 
















'w' 


X 


i— 


(XI 
















O'. 


II 


•X 


H- 


0 










1 — 1 




i — 


— ~ 


lU 


2T 
















X 


3 


Q 


liJ 


f— 










X. 








v~, 


Q 


< 














zr 


(XI 


#N 


»— 1 


_J 










< 




LO 


0 


vO 




U. 














X 


— « 


o 


vO 


(jC 










< 






LO 


o 


0 


CL 










0 






tO 


1 — 1 


0 


c-:) 










X 










rH 


K.J 










+ 






— 


h- 


h- 
















'1 


f— H 






r>' 




cr 










<T 


r 


t— < 




C 




f— 






In 




l— 


3 


(X 




!jL 




Ll- 




3 










Cl 


CL 






1 — 1 




<i 


11 




<i 


<T 






X 




X 




1 — 








1— 






^• 




u, 




X 


P 




X 


•< 










C 




< 


r) 




X 


X’ 






ll» 










lO 




* ■ 


X 




LO 


CL 




*' 




•' 






;;;7N 


!I 


o 


c:^ 




UJ 


r—* 




r-i 


(-^ 




u" 


— 


o 


c 






II 




II 




Li.. 


II 


3 


r— 1 


rH 


6 


V-. 


_J 


c 


\/ 


rH 


3 


— 








2 : 


h- 




• 




1 


2 : 




3 


K- 


}— 






vC 




r-- 


3 


f— . 




< 


2 : 


2 : 


< 


3 


r^t 


II 


r-^ 


+ 


H“ 


3 


2T 




♦ ( 


1 — 






> 




X 


Z' 


< 


6 


cr 


CL 


< 


< 


C 


X 


C‘ 


II 


0 


H- 


H-, 


Q. 


Cl 


c 


0 


C' 


u 


c 


3 


u 


DC 


IT) 



<x 

vO 

r 

(\i 

Cv 



c; 

C 

*> 

LiJ 



.0 

'''J 

\i 



iij ::; 



vO 

*■• 

r—i 

il 



o 

c 



U-J (X. 



Or' 

c 



C'J 

UJ 



r 

i( 

C 



<r 


3 


I— 




X 




1 — 






C' 


c 


• 






C'; 


•' 


i — 






*' 


0 








0 


* 




3 






Li 




L.j 




Ll _) 


3 




n 




3 


t— 


r'- 


' X. 




1— 


• — « 


0 


3 




►—1 


1— 


r~i 


•' 




1— 




3 


c 




•—4 


• 


<T 








3 


L J 


■'*< 




3 


Ll 


CL 


6 




L:.i 


CC' 


t-D 


z. 




C'' 


< 










3 


vO 


3 


(XI 


3 






< 


L3 


•> 


r 




3 


»— 1 


.0 




II 


L? 


Ln 




c. 




1—1 


X 


0 


0 




LT> 


l' 


0 






•' 


II 




•> 




tD 


3 




•, 


— 


1— 


Ll 


cv 




Ll 


U 


CC 




i— 


3 


f — 4 


< 


3 


DC 


1— 


2' 


3 


LO 


C' 


>— < 


(7) 






(r- 


1— 


,\/ 




<0 


1— 








3 


LL r— ' 




(/) 




LL 


v-H II 




i~ 


II 


I 


3 C 'i 


3 




r 


r 


t/ C' 








3 


\C X 


( ' 


T 




X 




1—1 






C' 


lO 






*.-> 


cn 


1— 


C 






3 


LL 


«■'{ 






ij_ 


V—* 


Ll 




h- 


f— . 


3 


CC 


<1 


cn 


X 


lT. 


«/> 


rr' 


< 


LO 


0 


i—H 




3 


i3 


Zl 






•' 




w 




3 


( — \ 




3 


LC 


3 




3 : 


<L 




<■ 


CO 


<f 


DC 1- 


<j" 


U 


€s 


QC 


(3 LT, ' 


II 




X> 


3 




3 


'.D 






3 3 


U' 




f'\j' 


3 


3 


cr. 


— 




3 


< i~i 


<< 


CXI 


rx 


< 


U 


3 






U 






Ljl' 

3 









(j 



* 



60 





























i — 




cr 




















1— 




1— 




\— 














c^. 




u. 


















H- 












c 




i — 




h~ 




H- 














a 






H 


C\. 






\L 








i- ■ 




u 












iLJ 




<iJ 




Ijj 














T~ 




«*T 


l 


c 


<x 




fy' 


(X} 


_J 


‘ . 


dc 


Cv.t 


'C> 






c\i 


1— 






cx 


cc 


r- 


Cl 


















'V' 




r -1 


_l 


~~L 


"C 


CO 


cc 


> 




(X) 


o 






•T> 




iL' 




ro 




c> 


cc 














• 






C.C 






1— 


<r 


T— 1 


c 


•X" 


< 


r— 1 


".■ 






r— < 


-J 






l~4 


< 


• 
















1- 


lt 


1'.. 


' r 


• 


o 




r.'i 




< — 1 


'Ll 


rn 




u_ 








loJ 


IJ.' 


•XJ 




o 


r- 
















• 


C' 


C\ 


ro 


r— 


o 


r'. 




K 








1— 


U-' 






h- 




1 — 
























Cl. 








r\' 


t-H 


i- 


pr 


Qr 


kC 


LiJ 


7C 


cv: 








fv; 




.J 


cz 


01 


























” 












f— 


•2T 








i> 






‘5 


>— 








Li 




r\j 


CM 


(\J 


fNJ 


r\i 


(N 


c\ 


c: 


( " 




(Xi 








fv 


a 


V— i 




Tn; 


•'f 


0 




Ov, 


u 


tx. 




(XI 




b 




c 


c 


c 


c 


c 


c 


c 


c 


O' 




c 


1 — 


K.f 




c 




c/:- 




c 


l.l 


or 




c 








c 








c 


c 


c 


c 


c 


c 


c 


I — 






c 


CL 






c 








c 




K 




c 


1.1 .- 


x> 




r 




c 




c 


, c 


c - 


c 


c. 


c.. 


c" 




O; 




c 


o 


V' 




c 


c::r 


C/") 






• 


„i 




c 


CO 


(<. 




c- 


.J 














1 — 1 


rH 




C-.0 


c. 




r* ^ 




r- 




r*H 




n. 








c 




r-H 








t-H 


X 






C- 


c 


c 


c 


c 


c. 


c 


l~ 


u 




r 


0'' 


c* 




c 




u 




cc 


b 






c 


o 


1 '*; 


a 


c . 




Ll.! 




c 


c. 


o 


C.' 


c . 


C' 


C- 


Cl 




_l 


c 


('O, 






c ' 








C ' 


I--I 






c 


-V- 


o. 




C •• 


b 


\~ 




c 


c. 


C ' 


r ' 


C- 


c 


c 


C' 


c 


Cl. 


( 




c 


_J 


c 




c. 






K 


X' 


o 


c 




X^' 




c 




1 — 1 




c 


c 


c 


c 


c 


c 


c 


c 


\C 




c 




if- 


n. 


c 


vX 


U ' 


V— . 


p- 




• 


llj 


c 


o 


c- 


-_J 


c 


'/) 


*7“ 




c 


o 


c; 


C' 


c 


c 


c. 


a 


c\: 




c 


o 


fx 




c 




OJ 


. J 


C' 


-J 


rc 


> 


c> 


lu 


(XI a 


o 




zc 




c 


c 


CJ 


c 


c 


r. 


c 


(f' 




OO 


c 




• 




o 


1— 




;■*> 


c 


t 




o 


r 


1— 






C' 


00 






c 


c 


c:. 


c 


c 


r 




• 


! 




c 


1~ 




LC 




<r 


• 




c 


o: 








1 — 1 


• 


'X 


r 


*■:* 






1 


rv.i 


f*-' 




ir\ 


vD 


1 — 


r*' 


(X. 


• 


c\: 


‘‘I 


“v 






..J 


rx 


X 


<1- 


o: 


<X' 


ll 


LO 




cx 


1/) 


o 


(*1 


fxl 


X' 


c 


c 


c 


(_ 


c 


c 


r. 


r-' 


a 


V-’ 


c 




cc 


CT; 


c 


L' 


^jD 


c/. 




o' 


cr 




c 




<0 




c. 




GO 


zz. 


c 


c 


C- 


c 


c 


(■ 


c 


c. 


r- 


i. ; ■ 


c 


.. 


r- 


u 


c 


rv 


r— < 






X 


.— 




.■r 


-J 


f- 4 




c 


> 




1— • 


c 


c 


c 


c 


c 


c 


c 


V.. 


• 


1 


c 




• 


ll 


c 






C / 


c 






z ^ 


c 






u- 


c* 






























Co 




l'^ 




c 


f- 






V. 


i — 


ll ■ 




o’* 


j— 


00 




r- 


h- 


r. 
















1. 


0 


fr 




c 


-V' 


o 






cc 








(k 


t' ! 




C.-X 


O' 


o 




c 


D. 






















o 




c.; 




(V 
















r"t . 




I’J 




(X 




• 






C\! 


c\ 


0*J 


r\i 


fV 


CXI 


(X 




cx 


1 — 


cx 




CX! 




(X! 


l: ‘ 


(X, 




(X 


(./ 


cx. 








CX! 




(XI 




(X 


(J) 


c 


c 


r 


c 


c 


c 


r 








c 




v' 




C” 


.X* 


X' 


KJ 




h- 




lC 






1 




c 




(. 




c* 


c 


o 


c 


c 


r 


c 








c 


f-- 


c 




r- 


''X 


.. ) 


I— 1 


c 




c 


’ < 


-• 


r>' 




•^r 


c 


• 


( 


u 


c 


c 


C.- 


c 


c. 


c 


c 


.. • 


: 


c 


c 


c 


r 




c 








c 


c; 




.z 


c 




C' 




( . 


01 


c 


UJ 


t—H 


r— : 


T— » 


1 — 1 


r— ' 




, — 


T'. 


r- 


c 


r-< 


c:> 


I — 1 




1 — ( 




r~ 1 




r-H 


> 


»-T 






oO 


f— H 


<r 


» — 1 


0^ 


r— 4 


X) 


c 


c 


c 


c 


c. 


c. 






r 


CX1 


c. 


H- 




o 


c 


L.J 


/•*' 


C' 


c 




1 


o'* 


V. 




c . 




c 


CO 


c 


o 


c 


c 


c: 


c 


c 


c. 


c 




c 




r 




c 


(\i 


C’ 


■•y 


o 


o 


c. 




c 


< 


< 




C' 




o 




cr 


cx: 


c 


c 




c 


c 


c 


c 


V ' 


c 


• 


’>. 


<f 


c 




C- 






iXi 


c 


1 


c 


1 < 


c 


o 


c 


• 


c 




r 


u 


c 


c 


C' 


c 


c 


c 


r 


O’', 


c 


f'L 


c 




c 


• 


c. 


<n 






r 




O- 




C ‘ 


(XJ 


< '' 


UJ 


c_ 


ll 


C‘ 


*— 4 


c 


c 


c 


c 


c 


c. 


c 


• 


c 


1 


c 


• 


c 


Q. 


c 


• 


o 


• 


c 


• 


r 




c 




X’ 


oO 


c. 


oo 


c 


' 


c 


c 


c 




c 


c o 




f— 


c 


Cl 


r 


H- 


c 


1 


c 


h- 


c 




c 




c 


r~~ 


c 


i — 




» — 1 


c 


♦— 4 


c 




c 


c 


c 


c 


c 


C r- 


c 


h- 


c 




c 


1— 




Cx, 


c 


f— 


V.C 


’l 


c. 


\- 


X 


i*“ 


c 


H* 


c 


o 


c 


Q 


c 




t-M 


r\j 


r' 


<r 




vTj 






1—' 






li! 


fX 




CO 


u- 


CO 


Q. 


<1' 


u 


*4" 


:.L. 


iC 


!i ’ 




zz 


vO 




o 


c 
















0' 




H- 








lC 




L>: 








C: 




0: 




o: 












(X 
















a 




i 




cc 












1 — 




6 




r\, 








.;J 




li.) 




















<1* 




c 




<r 








■cT 








•xC 








‘X 




f— 




H- 




• 
















T 




> 




cr 




o 




CL’i 




o 




Cf 




;'C, 




.0 




*— « 




»~4 




00 




















































*T* 
























































































CO 




“ 




Lf> 








vC 








pi 




r ■ 




r- 








rv' 



61 


















X 


























to 














(X 




X 






















< 














o 




o 


















• 














X 




X 




1— H 




CL 










u 




•vT 














lu 


X 


LiJ 




• 




CD 










O 




u 




o 










X 


t — t 


CD 




• 










r\j 


cfS 


< 




o 




f— . 


















• 








f-H 










X 
















X 


CD 




4\ 










Lu 


UJ 






iTi 




to 










h- 


CO 


X 












UJ 


X 


0- 


<. 


X 


\ 














< 




1— 1 




C\1 




w- 




CL 


< 


< 


1— 


X 


CO 




3 










X 


u 






■— 




X 




< 


H- 


I-* 


<1 


X 


X 




o 


• 








1— 


u 


LU 




X 




3 




f- 






O 


X 


X 




>- 


o 










•—1 


CL 




L3 




»— < 






X 


X 




<r 


Cl 






o 








o 
















X 


< 


< 


u_ 


X 






<L 


o 








to 




O 




*> 








< 


u 




c> 


u 






1— 


o 










o 


h- 








«— 1 




LJ 


»— 


►— * 






X) 






r-H 








LU 


UJ 






«— 1 








I 


to 


cO 


CO 


h- 






o 


CV 








CL 


DC 


o 








CL 




L9 


o 


O 


IL 


X 


u. 






ir\ 








< 


UJ 


X 




oc 




O 




O 


X 


X 


U 


XI 


o 




X 


i-H 










> 


*— ♦ 




CD 








u 






C 


X 






c 


«^ 










O 


> 




1— . 










cO 


to 


X 




r< 






r^ 








X 


u 


<■ 








X 




to 




(—1 


03 


U 1 


UJ 




to 


LTV 








b 


X 


X 




X 




lU 












~r“ 


D'i 




X' 


ti 










Cd 






X 




o 








— >• 


to 


H- 






u 


f\J 








< 




LU 




Q 




X 




— 


CD 




X 




Zj 




OJ 


\ 










LU 


CL 




X 




o 




o 


X 


o 




o 


X 






lU 






— 


<r 


03 


< 




o 








Uj 


o 


UJ 


>* 


X 






a 


\ 






O 


o 




H- 








X 




X 


X 


X 


X 


►— < 


to 












o 




C 






X 


* 


1— ( 




UJ 




U J 




X 






UJ 


C\i 






Cf'. 


X 


l__ 


1— 












o 


CL 


a 


1— 


X 






X 


\ 








o 




X 












X 


T*~ 


X 


X 


1 — ( 




• 


H- 


lO 






u 




cO 


xr 




1— 


X 


f— ( 




o 


lL\ 


o 


Ll i 


1— 


• 


>- 




f-H 






CJ 


to 




> 






•>— 






X 


1 — 




O 


X 




< 


X 


\ 






t— • 




u 


X 






X 


X 




3 








UJ 


to 


cc 


o 








»V 


X 


o 


CL 




X 


CD O 




— 




L.J 




c 




D 




\ 






— 


o 


u 


d... 




o 


1 — 1 








X 


0 


D 




o 


<r 


xi 


CO 






o 


u 


CO 








• 


CO 


«— * 


UJ 


X 


<]' 


Li ■ 




»— ♦ 




X 


\ 






o 


ro 




X 




tr» 


• 


\L 


X 


CL 


< 


d 


cr^ 


to 




LO O 


o 






<o 




X 


X 




LL 


• 


u 




< 


l~ 






K-. 


X 


X 




\ 








DC 


X 






U 




o 


X 


J— 




1— 


3 




o 


1— < 


o 


1— 1 






DC 


UJ 


1— 


X: 




o 


— 


_J 


CD 




X 


3 




X 




CD 










CD 


CD 








X 


(M 


03 




1— 


< 


X 




CD 


to 


X 


•— i 


*— • 






t — 1 


cc 


<r 


to 




X 






#> 


3 


X 


1— 


1— 


y—i 


X 


O X 


*- 








o 


y- 








X 


CO 


— 


X 


L 


3 


t—i 






-T- 










— 


UJ 


< 


X 




to 


CD 


< 


X 


X 


CO 


C 


iS> 


• 


< 












r— 1 


(Y 


CD 


1— 




< 




X 




1 — 1 








X 


X 


>- 




CD 






o 








• 


“T" 






X 










J J 




rc 


X 








c> 


O 


DC 


* 


V 




— 


UJ 














X 


<r 


X 


X 






r-H 


1— 


CO 


CL 


u 


LU 


r— i 


a 


1— t 


< 






I 


X 


o 


CL 




o 






•— 






UJ 


o 


a 




< 


•N 


X 






c. 


< 




o u 


CQ 




LU 


> 


CO 


>- 


o 


X 


<r 


X 


1— 










f- 




X 










U 


CL 


« — ( 


CTj 


X 


CD 


1— 


o 






o 






LJ 


0 


UJ 


SI 








ll 


;X 






o 








\— 




CL 






O 




CD 


0.1 


<c 






Ll. 


lu 




CD 




h- 


H“ 




ro 




a 








1— 


U" 


1— 












< 


< 


UJ 


X 


3 


> 


LU 










vD 


X 


3 




b 


H> 




X 




DC 


U.J 


> ll) 


0. 


no 


h- 




CO 






rH 


o 




<r 


o 


H> 






2: -o 


DC 




X 






3 


















X 


LU 






o o 




f— 






o 


O 




T- 






<T 


X 


UJ 


X 


X 


cc 




X 




(X 


to 


3 


e: 




LU 






i— 










X 


1 — t 








< 


uO 


Cl 


k—i 


u 


X 


LU 


■V: 


Lu 


— 








c/) 


Uj 


H“ 




UJ 


< 


• 


a' 


X 






Ui 


1— . 


X 


O 


u_ 


f\J 


-y 






I— . 


X 




>- X 


CO 


LU L9 


u» 


X. 


X 


c/} 


1— 


H- 


X 


H- 












i — 


cO 


X 


1— 




«v 


o 




1— . 


Cl 


X 






X 




X 








X 




1—. 


rr 






< 


CL' 


t— 1 


X 


< 


b 


Li- 




o 




CD 








X 


X 




X X 




«« 


CL 


o 


1— 


h- 


u 


LO 




X 




1— • 










o 


X 


sr 


H— ( 



^ -i 

a ^ 

o 

o f- 

U Li_ 



j — 




























(O 




1— ( 




























»V 
































m3 




cO 






to 






















r~H 




h- 


#x 




< 






















o 




3 


1— 




X 






















•> 




a 


X 
































Ui 




r— 






















o 




o 






t— 1 






















X 




X 


UJ 






















CO 




H- 




< 


1— 




h— 






















< 






< 




< 


















«N 




X 




to 


H- 




X 


















X 




3 




o 


cO 




i— 


















C<3 




’ r 




X 




























X 




LU 


XJ 




< 












X 






vO 




X 




rsi 


X 




h- 












a 






1 — 1 




)“ 






1— 




-r 












< 






o 








X 






3 










►-H 


1 — 














O 


> 
























-- 


3 






CD 




X 


• 










1 — 






H- 




X 










O 


to 








CO 


3 






< 




» 




UJ 


(D 






1— 








H-t 


X 














o 


X 




X 












y~ 






f-H 










(— ( 




u 


X 








to 


3 










— 




Xj 


3 




X 


X 








X 


o 




1 — » 


LU 




X 






X 














Ll 




— 


*N 


X 




3 




<r 


■> 




Q. 












X 






-t 


CO 


»— t 






•— 1 






< 








u 


o 




to 


1— 


►— t 






CO 






< 


X 








KW 












t.u 




XI 


X) 






3 








i — 


cO 






X 


X 


I— 




1— 






Q 


O 










CL 




vr 


o 


'O 


< 




< 






X 


X 








X- 


X 


n 


u 






X 




CL 








a 








r. 


o 




o 


Q 


vO 


o 




x; 


< 




X 










>—* 


d. 


LU 


X 


I.U 


r— 1 


X 






~r 


• 




UJ 










U- 




cc. 


CL 


o 






UJ 




to 


f — 


X 








c! 


X 


< 


\ 


X 




c 




o 




LU 


O 


h- 








o 




1— 


CO 


3 




H- 






to 




X 












iL. 


< 


UJ 


X 


X 






or 




X 


X 










CO 


O CD 


X 


u 




X 




<r 






< 










X 






Q. 


o 


< 


X 




X 


CO 


H" 


lJ 


f\j 








» 


X 


X 






u 






o 




■~y> 












<r 


UJ 


o 


< 


X 


o 


< 




o 


X 


b 


L- 


UJ 








a. 


r_. 


cc 


to 


o 


X 






X 


1_- 


UJ 


< 


X 










0: 


i— 




X 




I— 




X 








< 




— 




u_ 


o 




Ll 


fL 


o 


O 






• 


X 


3 


1— 




vO 




o 




6 


O 


UJ 


(— 


X 




XJ 




3 


3 




— 


t — ^ 






CD. 


UJ 












"T- 


Ll.' 


X 


;V 


X 


'■t 


o 




X 


U.J 




CL 


>- 


X 


nj 




1— 


X 


•V 


X 


o 








Ll. 


X 


H- 


Ll' 


1~ 


X 








< 








• X 


X 




X 


t—i 




X 


►—1 


X 


lU 




•X 


1 — 


1— 


:U 


X 


X 


0-'' 




;■ 


X; 


a 


r ~ 


X 


< 


Cl 








X 


X 


llJ 


CO 


• 




3 


X 




3 


< 


X 


< 




h-* 


h— 


3 


1— 


i — 




X) 






o 




:2: 


X 


3 


h- 




— ^ 


3 






H- 


CO 


r-* 




















X 


a. 


o 


»v 






c 




c.'- 


o 


O 


o 


o 


f— 1 


r-H 




I L 


L- 


X 


X 


X 






1—, 


X 


-I- 


X 


X 


X 


X 


X 




1 — 1 


3 


•— • 


lU 


Dc 


X 


X 


T 


U3 


C\( 


in 


CO 


vO 


o 


CD 






O 


X 


> 




CO 


0'' 


1— 


CO 




C\J 


nj 


ro 


r— 1 


m 




lU 






U.i 


>- 


— • 


— - 


— . 


— - 




— ' 




w. 


w 


■— 




X 


LU 


l-r 






1— 


1— 


K 




h“ 


i— 


1 — 


i— 


1 


1— 




< 


X 


o 


6 


< 


< 


< 


cr 


<L 


< 


< 


< 


< 


< 


< 




1— 


1— 


X 


X 


Ll J 


3l 




u 


jc. 






X 


MI 


rr 


7i. 


CO 






X 




lX 


X 


ft 


L> 


ai 


X 


X 


X 


DC 


ex: 


X 




X 


X 


< 


LL 


X 


CJ 




' — ^ 


o 


o 


o 


o 


CD 


o 


o 




o 


o 


u 




< 


X 


o . 


u 


X 


X 


X 


X 


X 


I.L 


X 


X 
































f— i 












1— 1 




o . 


<r- 


in 


CO 


cr> 


o 


r— 1 


C\J 





1—1 I — I r— I 



U u.u u u 



62 



O"! 

o 

t-H 

o 





i/) 




cn 


LiJ 




•—1 


•3 










X 


►-H 


u-^ 




h- 




r> 


■p' 


i- — 


vO 


6 




r-^ 




t— t 


C3 








X 


•X 




1) 





h— 


oc 




ft 
































C\J 


























< 


• 




































{XJ 


























































































0^1 


'vT 




#> 




) 




























C; 














Ijj 


















CT' 
































V- 














CC 












LJ 










_J 










































CY 












0.. 


0 




1— 




n-, 




























( 














0 




1 — \ 








<. 






X 




l~H 




























C j 


r**- 
















a 








i — 






> — t 


























'C* 








(X'. 


























U-1 




X 


















































3 








X 


h- 




Cl 


































0 
















3 








b 
































Cj 






— 


X 












U} 






















X 






















X 




















V 




0 








0 


0 






































0 












vJ 




X 








iij 


-J 




X 


























3 








uC 












x> 






r— 






Cj. 
























>- 








LO 






X 














X 




<r 








fv 


0 




<l- 


















0' 














O’ 


X 












X 




lo 








r) 


H~ 






















'Z 








0 






s 




7' 
















X 






u 






!— 




i — 














. : V 




<r 






•Cj* 




il) 


3 












X 




C'* 


rv 






u 


LU 




-y 




X 




c< 










X 




X 




X 


<T 




3 


•0 


fV 










LiJ 




b 


U 




vO 


0 


! 




»— • 




*-H 








Ov) 










C' 












— 


II 














t — 1 


iLi 




-t 








X 




X 




h~ 




II! 






— 




1 — 1 




:3 


C 






3 


X 














L/7 


U- ' 






0^ 


•< 




G. 




a 




S' 




G 






— 1 




rG 




~~} 


rr, 




r — 1 


— 












-i 




XJ 






T~? 


0 














— • 




<r 






— 
















(X 


vO 










X 




< 


r— 


0 


a■^ 


ul"' 


h~ 




*Sf 








or. 




L- 






X 




<1 




j 






> 


CC 












LJ 






r -1 


X 






0 












a 






c: 




.' c 




■J- 




L- 


cj- 






1 — * 










C\J 








•W^ 


<r» 


'0 


1.U 


X 
















X 


li 






-::t- 












• ' 




c 










X 




CO 


X 


O' 














— ^ 






X 


0 


fyy 




UJ 


»V 






r*' 






(J 


f— s 


3. 








0 


T“ 






rr 


<T 


— 


>- 


( — 1 




1— 














•-‘C 




:v“ 


PvJ 




— 




— 




7‘ _ 


1 — 1 












b 






X* 




r— 


f— 






SC. 




i 




e — 1 




, 


c. 


p. 






1 — 1 


JN 




f— . 


r. 






1* 








t— ( 






X 


Ll” 


CO 


I 1 




Ll? 




*—■ * 




C... 




-i 




> _; 












(<■' 








— 


t- 












rc 








<r 


s— 


OS 


CL 




o.-' 










C' 


X 




3 


— 




X 




LC 


>- 








-- 








LU 








— 




X 


<r 


< 




Q.' 




rx 




X 


<■ ‘ 


^ j 


* 1* » 


3 


1 — 1 


<t 


oz 


CC 


» — t 


0/ 








— ^ 








ul 


lO 




3 


■— 


X 


Cl 


n. 


i — 








— 








<’• 








C' 


z^~ 


cG 


X 




cr 


Ovi 




















•— 


<r 


zz 








i.A) 








it 


.•f 






(■ ) 


1 


0 


iu 




co 


hJ 


c 




CL 


o- 














0 






i7.- 


c> 


r~i 






T* 








1 — 


L J 




l~ 




<1- 


jZ 




L 


* ~ 


1 — 






r.; 












s 


X 


rv,! 




*r 


X 


X 




^ * 


• 


O’l 






.-f 




a\ 








1 


'~l_ 




1 


«> 




' -7 


.- • 








r\,- 


3 


«• 








*r 


vO 


0 




«x 


f~i 


C: 




(r, 


n 


, 


r-i 


CO 




t-H 




X 












i; 


!l 










r- ^ 


C\J 


X 


(M 




o-\ 


vn. 




*'L 


II 


•— 


Ll ■ 






X 


tl 




•— 




3 


3 




3 


vj: 


6 


r-l 




— 


r-l 


1 .J 


ijJ 


CC 


1) 


lG 


»— 1 


*» 


V 








•■V 


X 




ZD 


CiJ 


i 


e; 


3 






j__ 




3 


r— H 


— 


rH 






f-^ 




O'! 


5 


X 




V 






X 


w 


h- 


H- 




f~l 




r\ 


X 




c 


t— . 


' 


X 


X 


< — 1 


X 






X 




CC 


3 




V_- 




X 


1--1 


CM 




3 




»— t 


CC 


< 3 ; 


< 






0 






r— 


X 


X 


<3 


0 


LU 


-r?' 


0 


0 


X 


3 


I— 


L». ' 


X' 


>- 


X 


X> 


.-1 


LU 




<r\ 


0 


UJ 


;x 


0 




T* 




0 


r^) 


r. 


i — 


X 


h" 




CM 


h— t 


Ll. 


3 


t — t 


h- 




<— t 


X 


X 


1 


rs 


Cl 


1 — 


! — 




!U 


c\i 


t— t 


u 




*— » 




0^ 


h— 


< 




<r ' 




t— , 




»”- 




X; 


X 




— 










LU 


3 


«.u 


Z~ 




-V 


n 


Cl 




Lf) 


u 


w 




c 


6 




LU 


0 


L.^ 


0 


X 


b 


3 


0 


lO 


3 


LU 


LU 


V*-/ 


X 


i.U 


D-. 


■■ 1 


LiJ 


u" 


i.r 


'' ^ 


b 


?s 


•-L 


0 


X 





u_ 


X 


X 


u_ 


X 


DC 


0 


X 


u 


G. 


vu 


a 


Ci 


< 


a:> 


h-t 




LO 


X 




Q. 


iTj 


OC 


Z\ 


s 


C- 


X 


LiJ 


X 


0 


< 


X 


t — 1 


< — t 






1—1 




























































C<' 


<r 










C.‘ 












<r 


C'’ 


r-^ 


(\) 


CO 


C<^ 


c. 


,-- 


X 


X 


vO 




\j 










vO 


X 


a.i 


1 — \ 


f — ! 










X 












<r 


<L 




<T 


c<-^ 


<r 


C<^ 


nC' 




O' 






f\j 










vj- 


v^- 


vf 



U vj U u u 



63 



49 GO TO JUMP, (5?, 53) 

52 RFWIND 2 $ ASSIGN 53 TO JU:-)P S GO TO 45 

53 REWIND 2 i PRINT 12,IDR(K),< $ K=N S GO TO 23 























V.G 






















































cx 






















































o 














































o 


r— 1 






i— 














































vD 


<D 






o 














































o 


o 






3 














































r- 


i— 






•vf# 










on 




































o 


o 




















































o 


o 
















u : 








































!•: 




— 










3 








































rv 




> 






















<)- 






<!■ 






















f\I 




' — 










ll 












G 




(X 


G 




G 






> 


> 






t-H 








> 












c<' 










c 




c 


C 




C 








C’^ 






r:j_ 


'/ 


ll 














(3 












c 




c 


cu 




C 




r-^ 




> 






3 








/£ 






3 




-p" 


O 










r-H 






*-H 






vO 


lP, 


Li J 


u. 






3 




3 




UJ 


C 


c ■ 


3 




liJ 


3T 










c 




(X 


c- 




<x 






T'" 


X- 






3 


\./ 








v"’ 


3 


















o 




o 


O 




o 




o 




*' 








— 


rv 




«' 






Q 




trO 


3 










C- 




c 


c 




c 


v-f 












O 


3 


^0 






l: 


3 








'.J 










cx 




c\ 


rx 




cx 


vD 




r—' 


r-H 






i — 


O 








vO 


3 






C^i 


'3 










on 




0 


r' 




n' 




o 














O 






* 




o 
















c 


o 


c 


C' 


C 


C 


vO 


lo 


>- 


>• 






ro 




o 










LG 




h 












o 


h- 


r 




h— 


c 








Q- 


.r\ 




in 


•c^ 






fi-' 




O 






3 












<r 




<J 


iP. 




UG 




(-0 


I> 










r-i 












E'O 




i— , 


r\, 










c 


uT 


c 


c_ 


G* 


c- 






Uj 




OR 


’0> 










1 


— 








3 












c 


X’. 


c 


c 




e 


«~l 


V 


> 




I!" 


cr 




f— 










! 


NV 
















c: 


c 


c 


c 


C 


c. 




• 






r- 






■ y 


/ 




— - 












7L‘ 












b- 






c 




> 


— 


— 










— 


• — 




— 


>- 


>- 






X:' 
























o: 




t— i 




ir, 


.r. 


3 


ir: 


_J 




T— { 




Cl 


:s./ 




ll j 


3 


















«. 




1. 


— 


r 








<'!' 


G. 










Ci 








u. 












o 


<7 


<t 


• 


<j* 


■jj 


1 V 




i— I 


vC 








— « 






X 


;j 


3 






C'' 










x; 


o 


o 


c\> 


u 


fX 










1 ; 






f 






— 






n 
















r 


h- 


c 


c 


\ — 


c 


1 




r— 












li 


C' 


i 


1 


*1 


1 — 1 






> 










C 


L 


e 


c 


L'.. 


c 


— 




— > 














il 


r:) 




•— 
























r-. 






V 




o 






3 




' — 1 


— - 




o 


>/ 




o 


U 




r-H 




lO 








• 




c 


« 


fV' 


*— 


1 — 1 








3 






r— 1 


! 




w. 




r-i 


— > 


3 






^-4 




C- 


O' 




'C: 


o 


3 


c 






C'' 


Cr. 


1 — 




J 


3 


— 




cy. 


_i 


3 








GR 




Xf' 




CL 


c 


i— 


V. ' 


r 


1— 


C' 


n 




Li ■ 






O 




-3 


>- 


>- 


lLj 


c 




!— 




3 






4—4 




u. 


r-H 


Li. 


( — 


r-* 


U.' 


1— H 


1 — 1 




li_ 


L. 


3 


h- 


i~i 


1 — 1 


3 


O' 


ix* 






~~p 


h“ 




1—4 


li ■ 


!J- 




-T 


fX 




rv 


(X 




(X 


•— 






1 






~Z' 


iT 


:x 


X- 


Li_ 






^r-^ 




P,- 


3 


(3 




c 




L 




C' 


c 




e 


li. 


cy 


~ 




LL 


o 


L.J 


LU 


LU 


JJ 


3 


Li_ 




Zx. 


c 


3' 


LiJ 


3* 




J 




c 




c 


c 




c 




Cl 


rc 


cb 




eu 


on. 


3 




vx 


JTi 






Q 


3 


3 


q: 


LU 




j — 




r- 




r~H 


fX 




cx 








































3 


C 














o 


<r 


r-i 




vO 


CO 


<r 


i.n 


r~i 


C\j 


vD 


O 


r- 


Lf^ 












3 


3 














ir 


vO 


lO, 


vO 


lC 


eCR 




lG 


vO 


vG 


vO 


lO 


iG 


3 


cx 










1 i ; 






















































X 






















































Li.l 

















I 



64 



300 21 '^00 100 2 4 00 06''^ 0320001002 4 
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X ^cale-20 unitti ^inch. Y Seals- 2, unit/inch, 
figure 7- 253 cps so \are wave, 8 volts p-p. 200 usee, 
sampling. 
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X Scale-20 units/inch* Y Scale- 5 units/inch. 

Plgu3?6 6- Auto-e 0 in'’elatlon of Figure ? over 3900 samples. 
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X Scale- 20 imits/lnch. Y Soale-1 unit/lnoh. 
Pi^nire 9- 253 cps sawtooth wave, l+.o volts p-p. 

200 usee, sampling. 
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A Scale- 20 units/inch. Y Scale-1 unit/indh. 
Figure 10- Auto-correlatlon of Figure 9 over 3900 
aarr^ies. • 
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X Scale- 20 units/inch. Y Scale- 1 unit/inch. 
Figure 11- 253 cps sine wave at 6 volts p-p. 

200 usee, saij^llng* 
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X Scale- 20 u:iits/inch. Y Scale- % unit/lncli. 
Figure 12- Auto- correlation of Figure 11 • 

3960 sainplea . 
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Figure 13- Fluctuations about the mean of a D.C. v/ave* 
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X Scale- 20 vmits/inch. Y Scale- 1 unit/inch. 
Figure 15- Noise generator band-limited at 30 cps. 
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X Scale- 20 units/inch, Y Scale- 2 units/inch. 
i*'igure 17“ Noise genes»atop output. 2o6 usea» s&JtpiiiBg. 
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APPENDIX B 



NOISE REMOVAL IN CORRELATION, 

USING A PRIORI KNOWLEDGE OF NOISE 

Assume that as apriori information we are given -the expected form of 
the interfering noise, represented by a normalized correlation function 

formed in the absence of signal. We asjp|ime that when the re- 
ceived signal is corrupted by noise, the statistics of the received noise are 
described by this priori knowledge. Specification of the statistics of the 
noise by the correlation function implies that the noise distribution is ad- 
equately described by only it's first two moments. We shall assume that 
the observed date is of mean zero. 

Knowing the noise, it was shown that under cejjtain assumptions we 
could write 

^jL 

where %+N^^^ observed normalized coixelator output 

A 

Rss(J^ A T) is to be determined. 

A 

^NN^^^ T) is known apriori 

Obviously we must determine before Rss(k may be determined. 
Since the correlator produces a correlation function the author's first thoughts 
were to find 0 ( using the observed and known correlation functions. 

In the general case, the form of the correlation function of the signal 
is not know. For the class of input signals having a non-periodic correlation 
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function, any choice of which forces the output to be periodic, ie max 



R=R(0), would be fallacious. 

Therefore we shall only consider the class of signals having periodic 
correlation functions. 

♦* 

Let us consider the function 



F(^) = ](rS+n(0) - Max 

C Tjpo 






( 2 ) 



which may be written 






(3) 



F(o< ) = Max (l - *< ) - 
^ 4:0 I 

Since the functional forms of known explicitly 

no further analytic evaluation can be made of F(«^ ). 

Certain features of F(®^ ) are known however in the case of periodic cor- 
relation function. If FC*^ ) < 0, then ®< is too large,, and the output is not 
a valid correlation function since R(T4? O) > R(T = 0). If F(°<) > 0, not 
enough of the noise correlation has been subtracted, and the output will be 



excessively noisy for small T» 

The best choice of would then seem to be o( ^ such that F(o(^) = 0. 
The determination of o( must be done iteratively since no explicit 
analytic formulation ; for the f unction. j is known. 



I 

i 



78 



Consider Figure 1 below 




Evaluate F(^ ) and decide 

a) if F(<=^P > 0 

b) o< n = o< i if F(o^ i) < 0 j 

c) Terminate if [f(o^ pj < ^ 

Then 0<.^, = 

Having found such that - € £■ ^ ^ the portion of due 

to signal may be estimated by substitution in £q. (l)« that is 

Rss('^) = RiJn'Ct) - (4) 

At this point we have a combinational ratio which causes a periodic out- 
put. However there is no guarantee that the new correlation function formed 
in Eq. (4) is the result of a physically realizable linear system. This is be- 
cause the condition that the correlation be no larger for any T other than 
T = 0 is a necessary condition and not a sufficient condition. A sufficient 
condition for realizability is that its transform, the power spectral density, 
have no negative portions. 
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Rather than attempting to determine in the correlation domain, it is 
more fruitful to do a similar search in the frequency domain. The obvious 
advantage is that physical realizability is easily s^Jisfied. 

Transforming Eq. (1) we have 

PssOc ^ 

where the asterick indicates normalization. 

An iterative search, similar to the previous search, results in determina- 
tion of subject to physical realizability constraints. 

This is ensured by terminating the search when for any k, the value of 
Eq. (5) is less than .00001. 

If any value greater than this is used then the contribution will be neg- 
ative and realizability is violated. If a smaller value is used, more noise 
will be left than in the case 

This value of ©< = ^ ^ then used in Eq. (5) and in Eq. (1). 

The results of Eq. (5) were renormalized to have, unity area by the trap- 
ezoidal rule, that is 

L-1 

■^(kpi + kpj^) + 2 ^ kpi = 1 (6) 



from which we have 

k = 2/ (pj^ + Pl + 4 



L-1 

£ 



2 



Pi) 



(7) 



The results of Eq. (1) were renormalized by a factor 1/ R33(0). 

To evaluate the results assume that the correlation of the desired signal 
formed in the absence of noise Rs 3 (k^ T) is available for comparison. 
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An estimate of the noise power in the correlatc^T output may be found by 
evaluating 

No = W ^ - (1 - ^ c) T 

k=r ' 

The signal power in the correlator output may be evaluated several ways. 
One method is to assume that the scaled version of the pure signal, (1 - k ) • 

X 

Rss( 'I' = 0), represents the peak value of the periodic output. In the case 
of a single sin.usoid signal, the output power would be given by ((1 - ©< ) • 
RSS(0))V2. 

The input signal-to-noise power ratio may also be determined in several 
ways. One method is to compare the variance of the noise with that of signal 
plus noise to give 

(S/N)input = ^®S+N ■ 

An alternate method to determine S/N ratios and processing gains achieved 
is via the normalized power spectral density, obtained from the transform of 
the correlation function. This method is felt to be more accurate than deter- 
minations made from the correlation function because of the difficulties 
encountered in separating the signal components of. a correlation function 
from the noise components. 

Determination of signal to noise ratio from the normalized power spect- 
ral density proceeds as follows. 

The energy under the normalized density curve, with respect to a unit 
change of index i, is given by 
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X di = 1 



1 



Using the trapezoidal rule we may evaluate this as 

M 



^ ~ "7 ^M+1^ + 2 ^ Xj + 2X3 = 1 ; 

i=i2 

i^?S 



where Xg is the component due to the signal. 

The desired output, Y, is given in terms of frequency 
change of index and is given by 
Yi = kX^/E 
where k = 2 ^ tM 
= k'l 

A t = interval between samples in seconds 
M = the number of frequency estimates . 
Substitution of Eq. (1^ in Eq. (11) yields ? 

E = N + S 



where N = -^ (Xj^ +^M+1^ 2 



S = 2 E Y« 



M 

i=2 

i^S 



Xi 



The ratio of energy due to signal S, to that of noise, N is 

A 

= 2Ys 

1-2-i^ Af- 



S = 2 E^s 
N 



E - 2i Yc 



or in decibels 



(S/N)db = 10 log 10 (S/N) 



( 10 ) 



€ ( 11 ) 



rather than unit 

( 12 ) 



l4)3) 



(14) 



(15) 
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2. Comparison of Auto and Cross -Correlation 



No method of processing can improve the processing gain given by cross- 
correlation, which is taken to mean comparison of a signal plus noise with 
the completely specified input signal. This fact arises because of the a- 
priorii: information assumed to be available. 

When in general the desired signal is not known, auto-correlation, that 
is comparison of a waveshape with a time shifted replica of itself, may allow 
detection of the signal. In no case does it do it as efficiently as cross - 
correlation with an uncorrupted signal. 

If the noise can be assumed to be completely specified by it's correla- 
tion function, then it seems reasonable to assume that results approaching 
that of cross -correlation might be obtained by suitable removal of the known 
noise effects. 



Lee |_lj has developed models to estimate the processing gains possible 
for both auto and cross correlation, for the case of a siggle sinusoid 



For a single sinusoid input the desired output of the correlator ~ 

cos w, T where E = Em/ /~2^o 

At the input to the auto-correlator, let be th^ input noise in rms 
value, and let 



be the input noise-to-signal ratio. This definition while not conventional 
leads to compact formulation. 



2 

En\Sin(wt+ 9) corrupted by noise of zero mean and variance^T^ 




(16) 
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The output noise in rtiis' value may be shown to be 

1/2 • 



Noa = 






( 17 ) 



where N is the number of data samples considered. . For sinusoidal input, the 
rms signal output 

So^ = eV \pP ) (18) 

At the output of the correlator the signal to noise ratio is given by 



Roa = 20 1og^Q - 



oa 



(19) 



oa 






Substitution of Eq. (13), (14) and (15) into Eq. (16) yields 

If in the case of cross-correlation we assume t(ie signal plus noise is 
correlated with a local signal Em sin w t, the output noise in rms value . 

■' r . 

may be shown to be [_l]. 



N 



(|1 ^ e 2 S-J)] 






( 21 ) 



and the output signal rms 
S - ^ 

■ rv 



The gain in db, Rqc' is given by 



( 22 ) 



Rqc “ 20 log 



^oc 



1 0 Nqq 



db 



(23) 



This may be written using Eq. (16) 

N 



Rqc = 10 log 



1 + 25^1^ 

f' 

Inspection of Eq. (20) and Eq. (24) clearly shows the superiority of 
cross -correlation especially for largej’j^ (ie weak signals). 



(24) 
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( 25 ) 



The superiority is expressed by the difference 



G = V - "oa = 10 l°9l0 



• The concept of processing gain is perhaps more useful for comparison 
of any proposed method with these two standard methods. 

Let us define processing gain as the difference; in signal to noise ratio 



at the output and the input of the processor. 

f2 

= R . -in Inn. . ■ 

pa 



ie ~ 10 log^^Q ^2 



N 



7u 

2 



= 10 logio 



fi 



and 



'pc “ ^oc “ ^o^lO 
N f j2 



l + 4/i2 + 2j'i 

e2 



- 10 logio 



N 



4+ 



(26) 






db 



= 10 log 



1 + 2 






5 10 i°Oio f 



(27) 



Finally the difference in processing gains D, 



ac 



l^ac ~ ^pa ~ ^pa 10 l®9io (2 t^) db (28) 

The above forms are quite accurate forJ^j^ >3. However the quantity 

N. 

^runs; is not the most commonly used description of signal to noise 
ratio. Accordingly equations (26) and (27) have been recast in terms of in- 
put signal to noise power ratios in db, and plotted in Figure 2. Eq. (28) 
was similarly treated and is shown in Figure 3. 

Figure 3 conveniently shows the extent to which auto-correlation proces- 
sing is able to emulate cross -correlation processing. 
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BO- 





Figure 3*. RsocesstMG vs, SMR 
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Figures 2 and 3 are for 3750 samples < < 

3. Experimental Results 

The results obtained for the case of unity signal to noise ratio are ex- 
cellent, as may be seen by comparing Figures 5 and 10. Figure 5 is the 
power spectral density of a sine wave formed without noise. Figure 11 is 
the power spectral density of the 0 db signal to noise ratio input after the 
noise has been removed. A processing gain of 26 db was attained, as com- 
puted from the correlation function. This is mid-way between the performance 
of an auto-correlator and a cross -correlator. 

Since results were so encouraging, the -10 db signal to noise case was 
then considered. Several problems became evident that previously were mask- 
ed by the strong signal component. 

Determination of o{ becomes difficult at low signal-to-noise ratios due 
to the small signal component in the auto-correlation function, and led to 
the formation of physically unrealizeable correlation functions, as evidenced 
in Figures 13, 14 and 15. Inspection of Figure 12, the power spectral den- 
sity of the -10 db signal plus noise before noise removal, shows a strong 
component near 50 cps. It's transform, the correlation function (not shown) 
is difficult to interpret. 

This suggests that the signal-to-noise ratio measure of quality may not 
be the most meaningful for this type of signal detection. 

A more meaningful measure is the ability to distinguish a signal com- 
ponent in the power spectral density from a large noise spike. This ability 
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may be expressed in db as 



A = 10 log, « — for i s (29) 

^ max Yi 

If A is computed before and after noise removal, the increase in A is a 
direct measure of increased gain. 

In other words, at low signal -to-noise ratios the problem shifts from 
one of interpreting a correlation function, to that of evaluation of the psd. 

During the course of experiment it became evident that the spectrum of 
the noise may vary quite considerably from block to block. Since it was 
desired to remove the expected noise spectral density, more than a single 
block of noise must be used. By averaging the correlation function of noise 
over ten blocks much smoother spectra resulted (Figures 6, 7 and 8). The 
average noise correlation over ten blocks was used for the -10 db case. 

This of course implies stationarity over the period pf time represented by 10 
blocks of noise. 

The claim that averaging separately formed correlation functions is 
equivalent to an increase in integration time, when the process is ergodic, 
should be verified. Examination of Eq. (26) and Eq. (27) leads to the con- 
clusion that Gp^ should increase by 3db when the number of samples is 
doubled, ie 

^Gpa=10 logio ^ (30) 

where and N are the number of samples used. 

This goal was not reached, but was closely approached, falling short 

by about 1 to 1.5 db. 

The results obtained for the -10 db signal to noise ratio input are shown 
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in Figures 17 and 19, before noise removal, and in Figures 20 and 21 after 
noise removal. Noise removal had the effect equivalent to increasing the 
input signal to noise ratio by 14.46 db. Di scrimination against roise peaks 
was increased by 3. 42 db. 
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1074 FORMAT(45H IF IP50=1.PSJ DOME WHEN RJAU PLOTTED. IPSD= 13/) 

1075 FORMAT(45H IF IPRINT=1?°SD WILL BE PRINTED. IPRINT = 13/) 
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CALL TIMP 
PRINT 1012 
Ml=l % N2=6 
DO 700 J=l,21 
READ 701 , ( RSS( M) ,M = :-U ,.M2 ) 
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X Scale- l6 msec/inch. Y Scale- 0.5 units/inch. 
Figure Pure signal auto-correlation over 3750 

sarples. -i^stimated S/l\l ratio, plus 55 
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Figure 



Scale- 0.02 units/inch. Y Scale- 100 cps/inch, 
5- Pure signal power spectral density over 
37^0 saraples. 125 estimates 5 cps apart. 




X Scale- .001 units/inch. Y Scale- 100 cps/inch. 
Figure 6- Power spectral density of noise, 3750 ssjnples, 
125 estimates 5 cps apart. 




X Scale- .001 -units/inch, Y Scale- 100 cps/inch 
Figure ?- ?ov/er spectrs.1 density of noise, l8?p0 samples, 
12^ estimates 5 cps apart. 




X Scale- »001 units/inch. Y Scale- 100 cps/inch. 
Figure 8- Power spectral density of noise, 37^00 samples, 
125 estimates 5 cps apart. 
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X Scale- li.O msec./iiich. Y Scale- 0.5 uuits/inch. 
Figure 9 - Auto-correlation of ODB. signal plus noise. 
3750 samples . 
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I 




X Scale- 0.02 units/inish, Y Scale- 100 cps/inch. 
Figure 10- Power spectral density of ODB. signal after 
noise removal. Alpha calculated on RBAR(T). 
3750 samples . 
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-323^ "01S -019 -005 009 _ . 013 




X-Scale- ii-0 msec, /inch. Y Scale- ,5 units/inch* 

Figure 11- ODB. signal after noise removed. Alpha 
determined on RBAR(T), 3750 samples . 
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X Scale- 
Figuee 12- 



^ Scale- 100 cps/inch. 
~ 10 uo. psd. before noise removal. 

125 estimates, 5 cps apart, over 3?50 samples. 
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X Scale- .0J>0 units/inch. Y Scale- 100 cps/inch. 
Figure 13- -lODB. psd. after noise removed* 37^0 samples* 
Alpha determined on RBAR(T). 
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PiSu^o^Ju-® noLe^remOTaJ? !ll^o"s^lea 

Alpha determined on RBAR(T). ^ ample s. 
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. X Scale- .02 unita/inch, ■ Y Scale- 100 cps/inch. 

Figure 15- -lODB. psd. after noise removal. 18750 samples* 
Alpha determined on RBAR(T). 




’ X Scale- ii-0 msec. /inch. Y Scale- 0.5 units/inch. 

Figure l6- -lODB. signal after noise removal. 10750 sample 
Alpha determined on RBAR(T). 
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X Scale- 16 msec/inch. Y Scale- 0.5 units/inch, 
J’igure 17- -lODB. auto-correlation before noise 
removal. 22500 samples. 
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X Scale- 16 msec/lnch, Y Scale- 0.5 units/inch. 
Figure 18- Estimated auto-correlation of noise in -lODB. 
signal plus noise input, averaged over 
37500 samples. 
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X Scale- .002 units/inch. Y Scale- 100 cps/inch. 
Figure 19- -lODB. power spectral density before 

noise removal, 22500 samples, 125 estimates, 
5 ops apart. 
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X Scale = 16 insec/inch. Y Scale - .5units/inch. 
Figure 20 - -10 DB. Signal afcer nois® removal in 

frequency domain. 22500 samples. 
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X Scale = .02 units/inch. Y Scale - 100 cps/inch. 
Figure 21 - Power spectral density of -10 DB. SNR 
after noise removal in frequency domain. 
22500 samples. 
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APPENDIX C 



DETECTION AND LOCATION OF SIGNALS IN NOISE 

In Appendix A and B methods were developed for digital analysis of 
analogue signals. 

The data utilized in the above appendices was in a sense artificial, in 
as much as it was generated in the laboratory. » 

To test the proceedures developed using a live signal, where signal 
characteristics and noise statistics are generally not known, a 200 cps 
sinusoidal signal source was placed in an unknown noise field produced by 
patrons, of a cafeteria. 

The signal source was set up at one end of the, dining room, and its' 
volumn raised until complaints were voiced. The source level was then 
lowered until the complaints ceased. This level, which we shall refer to as 
the 0 db level was measured to be 38 milliwatts at. the voice coil of the loud- 
speaker. At this level the sound was occasionally audible at various loca- 
tions within the room, but was not detectable aurally at the microphones. 

A stereo recorder was set up at the far end of the room, using two micro- 
phones symmetrically placed about the longtitudinal center line of the room. 

A number of three minute recordings were made of signal plus noise, re- 
duced signal plus noise (-9.3 db), and noise alone,. 

The noise level of the crowd varied very considerably over the hour and 
was punctuated by an occasional dropping plate. The period from 12:30 to 
12:40 was selected for further study as the number of diners was essentially 
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constant over this period. The details are listed in Table L 



By cross -correlating the left and right channel information it is possible 
to locate the source and to estimate its frequency. . 

The bearing of a source, with respect to the center-line may be deter- 
mined by measurement of the difference in arrival times of sound waves at 
the receiving sensors. 

Assume as in Figure 1, a distant sinusoidal source, at^amangle 0 mea- 
sured from the center-line or dead-ahead position, and two sensors symmet- 
rically placed about the center-line spaced d units spart. 

Assuming plane wave propagation, with no multipath effects, and in the 
absence of noise, the difference in travel time or phase of ray 2 is a direct 
measure of the angle 0. 



Figure 1. (a) The basic DF Problem, (b) Cross -Correlation of M]^, M 2 . 





.T/2 



Forming Rj^^ , lim T-^<=Q 



( 1 ) 
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/W 

A [cos f + ( y) dt 




o 



cos[Wpr - arjdt (2) 



= 0 + cos (WqT -y) 



(3) 



The first peak in the correlation curve occurs when 




(4) 



Subsequent peaks in RjvilM2 occur when 

WoO* - y = N (2TT ) ^ N = 0, 1, J2. 3 ... 

Peaks obtained for n H 0 represent ambiguities in the bearing angle. The 



The first peak at ^ is of interest. In this case using room measure- 
ments and an assumed velocity of propogation of 1129 ft/second, % was 
calculated to be 1. 5506 msec, for an angle of 10. 4°. For small angles we 
may assume the time and the associated angle are linearly related. Under 
this assumption the arrival time would change very nearly 150 jX: sec. per 
degree. 

The available processing system however could sample the record at 
best every 200 :p. sec. , resulting in a bearing accuracy of about 1. 33®. To 
improve the bearing accuracy, the data was slowed by a factor of four, 
grequency shifting all components down by a factor of four. This has the 



first ambiguity occurs when Tp = — sin ^0orA0 = 0^-0 = sin”^ 
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effect of increasing the arrival time by a factor of four resulting in a bearing 
accuracy of 1. 33/4° or . 333°. 

The mechanics of the speed reduction are shown in Table 2. 

This is more bearing accuracy than is justifiable but since a reduction 
in data speed of at least two was necessary to digiti 2 Bja 5 Kc analogue re- 
cord, the extra factor of two was used to increase bearing accuracy. 

The data speed reduction meant that the original. 200 cps signal would 
appear as 50 cps with the first peak of the cross -correlation curve at 'T = 

6. 02 msec. The data was sampled at intervals of 205 jLLsec. , and hence 
the signal, if present, should occur between 29 and 30 shifts of the correla- 
tion function argument. 

2 . Data Analysis 

Twenty blocks of data were digitized for each of three cases; signal to 
noise ratio 0 db; -10 db; and noise alone. 

Due to difficulties experienced with the clocking system, only the first 
block of a twenty block run was clocked, with clocking of subsequent blocks 
internally. Due to the variations in the time necessary for the 163 tape uhitL 
to recover from a writing operation, a cumulative error in clocking arose: 

Hence only the first few blocks are time synchronized to an acceptable degree. 
Fortunately this was long enough for sucessful signal detection. 

The data was normalized in the computer by the computed estimator of the 
standard deveation, and the cross -correlation of the normalized data formed. 

The time of integration was fixed due to the finite block length. Assum- 
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ing a process is stationary and ergodic, (at least short-time stationary) the 
effective Integration time may be extended by averaging in time, ie 



N 



^ ^ N=1 



(5) 



- K 

and Riab(T) = — Riab(T) 

K i=i 




i=l N=1 




( 6 ) 



The running average cross -correlation function of the 0 db signal plus 
noise was formed over the first five blocks and is shown in Figures 3 through 
5 as a "hatched-in" curve. 

3. Results 

Experiments with averaging more than five curves were not successful, 
due to the previously mentioned cumulative timing error. 

The cross -correlation of 0 db signal plus noise was abserved to be es- 
sentially sinusoidal with the first peak occurlng between 29 and 30 shifts 
exactly as predicted. 

The frequency of the sinusoid however was observed to be 150 cps, or 
three times the expected signal frequency, indicating detection was being 
made on the third harmonic of the source! 

No firm explanation for this is offered. However it is noted that the 
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oscillator produced 0. 1% third hormonic, and the 20 push-pull amplifier 
used would be expected to produce considerable cross-over distorion when 
operated at such low power levels. In addition the tape re-recording process 
tended to favour the upper registers. 

A power spectral density was made of the 0 db signal plus noise, and 
of noise alone, using a 6. 125 cps filter bandwidth.. The normalized results, 
shown in Figure 2, reveal the pronounced 150 cps component in the signal 
recording, with the general shape of the remaining portion being similar to 
the noise power spectral density. 

The programming mecessary is included as Program SIMSIG in Appendix 
B. 

The results obtained with the signal at -9.3 db were inconclusive and 
are not included. 

4. Conclusion 

Provided accurate time synchronization between right and left channels 
can be maintained during the digitizing process, cross -correlation of the two 
channels provides both bearing and frequency information. 

The process of averaging cross -correlation curyes must be done with 
care, ensuring accurate time synchronization and continuing the averaging 
process only as long as the signal characteristics remain stable. 
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ANALOGUE E)ATA DIGITAL 
TIME STATE RUN 



I2;30-12:33 "O" DB. 


F 


0001 


0021 0001 


0024 


(R) 








TO 










0001 


0021 0024 


0024 








0001 


0032 0001 


0024 


(L) 








TO 










0001 


0032 0024 


0024 




12:35-12:33 "-10"DB. 


J 


0002 


0021 0001 


0024 


(R) 








TO 










0002 


0021 0024 


0024 








0002 


0032 0001 


0024 


(L) 








TO 










0002 


0032 0024 


0024 




12:41-12:44 NOISE 


C 


0003 


0021 0001 


0024 


(R) 








TO 






Digitized on Reserve Tape 18 


0003 


0021 0024 


0024 




Sampling Interval 205 ^s. 














0003 


0032 0001 


0024 


(L) 








TO 










0003 


0032 0024 


0024 




TABLE I: Details of Data Digitized. 









STEP 


PLAYBACK 


BAND-W. 


RECORD 


BAND-W. 




Concertone 
7.5 IPS 


50 cps 
10 Kc. 


Ampex FR-100 


0 - 1. 1 Kc 




Ampex 
7.5 IPS 


0 - 650 cps 


Concertone 
7.5 IPS 


20 - 650 
cps 




Concertone 
3.75 IPS 


10 - 325 cps 


Ampex 
7. 5 IPS 


10 - 325 
cps 



TABLE II; Frequency Reduction Process. 



SIGNAL 

200 

cps 

100 

cps 

50 

cps 
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X Scale- 2.05 msec/inch. Y Scale- 0.5 units/inch. 
Figure 3- ODB. signal cross-correlation over 3950 
sanples. 
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X Scale- 2»05 rasec/inch. Y Scale- 0,5 
Figure I;.- 0D3. signal cross-correlation over llo^u 

•• samples. 
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X Scale- 2.05 msec/inch. Y Scale- 0.5 units/inch. 
Figure 5- ODB. signal cross-correlation over 19750 
samples. 
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X' Scale- 10.25 msec/inch. Y Scale- 0.1 volts/inch, 
Figure 6- ODB. signal aut'o-correlation over llo^O 
samples . 
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^ Scale- 10,25 msec/inch. Y Scale- 0,5 volts/inch. 
Figure 7 - Noise auto-correlation over ll 850 samples. 
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-’igure 6- Relative pov/er spectral densities of noise 
alone, and signal .plus noise. Normalized 
to 100/» . 
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